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CHAPTER  I 
BASIC  CONCEPTS 

The  concept  of  inertial  navigation  is  not  new.   Sir  Isaac 
Newton,  in  his  "Mathematical  Principles  of  Natural  Philosophy," 
gave  the  laws  of  mechanics  on  which  self-contained  inertial  navi- 
gation is  based.   At  the  time  of  the  publication  of  Newton's 
treatise  there  were  no  means  available  for  sensing  accelerations 
accurately  or  for  economically  determining  positions  if  sensing 
had  been  possible.   It  has  taken  almost  three  hundred  years  to 
develop  the  necessary  hardware  and  software  to  make  self-contained 
inertial  navigation  a  reality. 

A  brief  description  of  an  inertial  navigation  system  is  that 
we  have  a  "black  box"  which  senses  both  accelerations  and  mass 
attraction.  We  eliminate  the  effect  of  the  mass  attraction  and 
then  use  our  initial  position  and  velocity,  elapsed  time,  and  the 
double  integral  of  the  acceleration  to  determine  an  updated 
position.   In  practice  this  is  done  many  times  a  second,  provid- 
ing a  continuous  readout  of  incrementally  corrected  position. 
Henry  described  part  of  the  mechanism  in  the  following  words: 

A  basic  requirement  of  an  inertial  navigation  system 
is  that  the  acceleration  vector  be  known  at  all  times. 
Knowledge  of  the  direction  of  this  vector  requires 
either  maintaining  an  acceleration  sensitive  element 
in  a  fixed  orientation  or  continuous  knowledge  of  the 
orientation  assumed  by  the  acceleration  sensitive 
element. 


The  system  may  take  any  one  of  the  various  forms  in 
which  the  acceleration  sensitive  element  or  elements 
are  maintained  relative  to  the  local  gravity  vector 
(geodetic  level)  or  geocentric  level  (that  is,  pointed 
towards  the  center  of  mass  of  the  earth).   The  stable 
element  may  be  fixed  with  respect  to  the  earth,  ro- 
tating around  the  earth's  axis  at  earth  rate,  or  may 
simply  be  fixed  in  inertial  space.   An  azimuth  line  on 
a  level  platform  may  be  oriented  to  point  at  geodetic 
north  or  in  some  other  arbitrary  direction,  (k) 

During  the  1950' s  and  1960's  a  tremendous  effort  has  been 
devoted  to  eliminating  errors  in  sensing  systems  and  servo  sys- 
tems so  that  the  acceleration  sensitive  elements  may  be  maintained 
in  a  known  orientation.   As  stated  by  Moore: 

The  development  of  inertial  navigators  has  been  one  of 
the  major  technological  achievements  of  the  twentieth 
century.  The  outward  simplicity  of  the  mechanization 
of  these  navigation  systems  and  the  ease  of  statement 
of  the  simple,  fundamental  scientific  principles  on 
which  they  are  based  mask  a  major  scientific  engineer- 
ing and  manufacturing  development  program.   It  has 
taken  the  combined  efforts  of  engineering  specialists 
in  electronics,  mechanical,  chemical,  and  electrical 
engineering  as  well  as  physicists,  mathematicians, 
metallurgists,  metrologists,  skilled  machinists,  and 
experienced  electronics  and  electromechanical  manufac- 
turing departments  working  in  close  coordination  to 
develop  and  produce  the  apparently  simple  end  product.  (6) 

It  has  become  apparent,  in  searching  the  literature,  that 
some  degree  of  confusion  exists  about  our  ability  to  eliminate  the 
effect  of  mass  attraction  from  the  values  sensed  by  the  "black 
box."   Statements  such  as  that  of  Henry,  quoted  above,  equating 
the  gravity  vector  with  the  normal  to  geodetic  level  are  common. 
In  those  cases  where  some  acknowledgment  is  made  that  the  magni- 
tude and  direction  of  gravity  do  not  agree  with  simple  formulas, 
the  effect  is  considered  as  random  noise  in  the  system. 

It  is  the  purpose  of  this  paper  to  assess  the  impact  of  the 


3 
variations  of  gravity  on  a  vessel,  such  as  a  survey  ship,  which  is 
likely  to  remain  within  a  limited  area  for  an  extended  period  of 
time,  such  that  the  effects  of  anomalous  gravity  are  systematic 
rather  than  random,  and  to  suggest  some  system  to  eliminate  these 
effects  from  the  navigation  solution.. 

The  effects  of  gravity  on  the  three  forms  of  the  system  as 
described  above  by  Henry  will  be  investigated  to  see  if  the  dif- 
ferent forms  are  affected  differently.   These  three  forms,  in 
order  of  increasing  complexity,  are; 

a.  A  system  in  which  the  orientation  of  the  sensor  axes  is 
fixed  in  inert ial  space, 

b.  A  system  in  which  the  orientation  of  the  sensor  axes 
maintains  a  fixed  relationship  to  the  earth,  with  the 
acceleration  sensitive  elements  rotating  at  earth  rate 
about  an  axis  which  is  parallel  to  the  earth's  spin  axis, 
and 

Co   A  system  in  which  the  orientation  of  the  sensor  axes  is 
dependent  on  the  navigator's  position,  where  one  axis  is 
aligned  with  the  normal  to  the  adopted  ellipsoid  at  that 
position,  one  axis  points  north  in  the  plane  tangent  to 
the  ellipsoid,  and  the  third  axis  completes  an  orthogonal 
system. 
Three  coordinate  systems  will  be  employed  in  this  study,  all 
of  them  being  right-handed  orthogonal  systems.   The  X  system  is  an 
earth-fixed  system,  with  the  X3  axis  parallel  to  the  earth's  ro- 
tational axis,  positive  to  the  north,  the  X°  axis  parallel  to  the 


Ya 
Xa 


Y0' 


Figure  1. 
Relationships  between  coordinate  systems 


equatorial  plane  and  the  Greenwich  meridional  plane,  positive  in 
the  direction  from  the  center  of  the  earth  to  the  zero  meridian, 
and  the  X1  axis  positive  in  the  direction  from  the  center  of  the 
earth  to  the  ninety  degree  east  meridian. 

The  Y  system  is  fixed  in  inertial  space  with  respect  to 
orientation,  with  Y2  =  Xa ,  and  is  related  to  the  X  system  such 
that,  at  time  t  =  0,  Y°  =  X°  and  Y1  =  X1  ,  and  at  any  other  time 
there  exists  an  angle  of  fit  between  the  respective  Y  and  X  axes  in 
the  sense  shown  in  figure  1. 

The  Z  system  is  a  function  of  position,  and  is  oriented  to 
the  ellipsoidal  normal  such  that  the  Z°  axis  is  parallel  to  the 
normal  and  is  positive  in  the  direction  of  the  outward  normal 
(generally  up),  the  Z1  axis  is  positive  in  the  geodetic  east  di- 
rection, and  the  Za  axis  is  positive  toward  geodetic  north.   In 
figure  1,  the  Z  system  is  shown  as  being  topocentric  while  the  X 
and  Y  systems  are  shown  as  being  geocentric.   This  distinction  was 
made  in  the  interest  of  clarity  of  the  figure  only.   All  three 
systems  will  be  used  as  geocentric  systems  so  that  no  translation 
between  origins  is  required. 

Rotation  matrices  used  to  rotate  the  coordinate  systems  about 
any  coordinate  axis  will  be  of  the  form  R.(q)  where  the  index  "j" 
indicates  the  axis  about  which  the  rotation  "q"  takes  place,  with 
the  sense  of  "q"  in  accordance  with  the  right-hand  rule.   Thus 
^  =  Ra(-0t)X  and  X  =  Ra  (-\  )Ri  (<p)Z  where  the  position  vectors  are 
defined  as  columns. 


The  following  symbols  will  be  used  throughout  this  study: 

a  -  semimajor  axis  of  the  ellipsoid 

b  -  semiminor  axis  of  the  ellipsoid 

e  -  eccentricity  of  the  meridian  ellipse  =  (a3  -  ba )  /a 

— ♦ 

g  -  gravity  vector 

h  -  ellipsoidal  height 

1  -  unit  vector  along  the  zero  axis  of  a  coordinate  system 

j  -  unit  vector  along  the  one  axis  of  a  coordinate  system 

a* 

k  -  unit  vector  along  the  two  axis  of  a  coordinate  system 

k  -  Newtonian  constant  of  gravitation 

G  -  gravitation  vector 

H  -  orthometric  height 

M  -  radius  of  curvature  in  the  meridian  =  a(l-ea )/tia 

If  -  radius  of  curvature  in  the  prime  vertical  =  a/W 

as  a  subscript,  N,  indicates  navigator 

R  -  radius  vector 

W  -  (1  -  essin2tp)^ 

or  -  azimuth 

V  -  value  of  normal  gravity 

Q  -  geoid  undulation 

tj  -  prime  vertical  component  of  deflection  of  the  vertical 

\  -  geodetic  longitude,  positive  east,  measured  from  X°  axis 

\'  -  astronomic  longitude,  positive  east 

£  -  meridional  component  of  deflection  of  the  vertical 

tp  -  geodetic  latitude,  positive  north 

<p  -  geocentric  latitude,  positive  north 
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cp'   -  astronomic  latitude,  positive  north 

u)   -  generalized  rotation  vector 

A   -  inertial  longitude,  positive  east,  measured  from  the  Y° 

axis  =  \   +   fit 
Ci        -      earth  angular  velocity  vector 
()   -   first  derivative  of  ()  with  respect  to  time 
()   -  second  derivative  of  ()  with  respect  to  time 
fig  -  gravity  disturbance 
Ag  -  gravity  anomaly 

0°  -  unit  vector  in  the  direction  of  ( ) 

Vector  quantities,  other  than  i,  j ,  k,  when  written  without  the 
sign  of  the  vector,   ,  denote  magnitude  of  the  quantity  considered. 

There  are  several  quantities  used  in  this  study  which  require 
specific  definition  in  order  to  avoid  a  chance  of  confusion. 

Meridian  plane  -  The  plane  containing  the  normal  to  the 
ellipsoid  at  a  given  point  and  the  axis  of  rotation  of  the  ellips- 
oid.  The  meridian  is  the  intersection  of  the  meridian  plane  and 
the  ellipsoidal  surface  for  that  half  of  the  ellipse  which  con- 
tains the  given  point,  and  is  bounded  by  the  intersection  of  the 
ellipse  with  the  rotational  axis.   Meridian  and  meridian  plane  are 
frequently  used  interchangeably. 

Prime  vertical  plane  -  The  plane  containing  the  normal  to  the 
ellipsoid  at  a  given  point  which  is  perpendicular  to  the  meridian 
plane  at  that  point.   The  prime  vertical  is  the  intersection  of 
the  prime  vertical  plane  and  the  ellipsoidal  surface.   It  is  tan- 
gent to  the  parallel  of  latitude  at  the  given  point.  Prime 
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vertical  and  prime  vertical  plane  are  often  used  interchangeably. 

Latitude,  geodetic  -  «p  -  the  angle  between  the  normal  to  the 
ellipsoid  at  a  given  point  and  a  plane  perpendicular  to  the  axis 
of  rotation  of  the  ellipsoid. 

Latitude,  astronomic  -  cp'  -  The  angle  between  the  gravity 
vector  at  a  given  point  and  a  plane  perpendicular  to  the  mean 
rotational  axis  of  the  earth. 

Latitude,  navigator's  -  <p»  -  The  value  which  the  inertia! 
navigator  indicates  as  being  descriptive  of  one  component  of  its 
position. 

Latitude,  geocentric  -  <p  -  The  angle  between  the  radius 
vector  to  a  given  point  and  a  plane  perpendicular  to  the  axis  of 
rotation  of  the  ellipsoid. 

North  error  -  The  difference  between  either  the  geodetic  or 
astronomic  latitude  and  the  navigator's  latitude,  converted  to  ay- 
linear  measurement.  Whether  the  error  is  computed  for  the  geodetic 
or  astronomic  latitude  will  be  indicated  in  the  text. 
North  error  =  (cpM  -  <p)M,  or  North  error  =  (<pM  -  <p'  )M 

Longitude,  geodetic  -  V  -  The  angle  between  the  meridian 
plane  at  a  given  point  and  the  meridian  plane  of  Greenwich. 

Longitude ,  astronomic  -  X '  -  The  angle  between  the  astronomic 
meridian  plane  of  a  given  point,  and  the  astronomic  meridian 
plane  of  Greenwich,  where  the  astronomic  meridian  plane  is  defined 
as  that  plane  containing  the  gravity  vector  at  a  given  point  which 
is  parallel  to  the  axis  of  rotation  of  the  earth. 

Longitude,  navigator's  -  \H    -  The  value  which  the  inertial 


navigator  indicates  as  being  descriptive  of  a  second  component  of 
its  position,  in  a  direction  orthogonal  to  the  navigator's  lati- 
tude o 

East  error  -  The  difference  between  either  the  geodetic  or 
astronomic  longitude  and  the  navigator's  longitude,  converted  to 
a  linear  measurement ■  Whether  the  error  is  computed  for  the  geo- 
detic or  astronomic  longitude  will  be  indicated  in  the  text. 
East  error  =  (,\N  -  \)Ncostp,  or  East  error  =  (\N  -  \')Ncos<p 

Deflection  of  the  vertical  -  £,,  r\   -  The  angular  distance 
between  the  astronomic  position,  q>',  \',  of  a  point,  and  the  geo- 
detic position,  (p,  \.,  of  the  same  physical  point.   The  deflection 
is  divided  into  a  component  in  the  meridian  plane,  £,,  aQd  a  com- 
ponent in  the  prime  vertical  plane,  i\* 

£  =  (p«  -  <p 

H .  ■  (X*  --   \)coscp 

Local  vertical  rate  -  The  time  rate  of  change  of  the  local 
vertical..   This  is  divided  into  components  in  the  meridian  plane 
and  in  the  prime  vertical  plane.   For  an  ellipsoid,  the  meridional 
local  vertical  rate  is  the  rate  of  change  of  the  ellipsoidal  nor- 

o 

mal  in  the  meridian,  or  the  rate  of  change  of  latitude,  <p.   For 
the  geoid,  we  must  consider  the  latitude  rate,  plus  the  rate  of 

change  of  the  meridional  component  of  the  deflection  of  the  ver- 

.    »    . 

tical,  or  <p  +  £  =  <p'  » 

Similarly,  the  prime  vertical  local  vertical  rate  for  the 
ellipsoid  is  \.coscp,  while  for  the  geoid,  the  rate  of  change  of  the 
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deflection  of  the  vertical  component  must  be  considered,  and  the 
prime  vertical  local  vertical  rate  is  \cos<p  +  r\   =  \'cos<p.  When 
this  term  is  used  in  the  text  it  will  always  refer  to  the  geoidal 
application  rather  than  the  ellipsoidal.   The  meridional  local 
vertical  rate  is  also  referred  to  as  the  astronomic  latitude  rate, 
while  the  prime  vertical  local  vertical  rate  divided  by  the  cosine 
of  the  latitude  is  the  astronomic  longitude  rate. 

Gravitation  vector  -  G  -  The  force  exerted  on  a  body  of  unit 
mass  because  of  mass  attractions   For  a  central  force  field, 

G  =  -  -j3-  ,  where  m  is  the  mass  of  the  attracting  body. 

Gravity  vector  -  g  -  The  combined  force  of  gravitation  and 
the  centrifugal  force  caused  by  the  rotation  of  the  earth,  when 
applied  to  a  body  of  unit  mass. 
g  =  (!  -  f5x(?Jxl) 

Height,  ellipsoidal  -  h  -  The  linear  distance,  measured  along 
a  normal  to  the  ellipsoid,  from  the  ellipsoid  to  a  given  point. 
Height  is  positive  along  the  outward  normal o 

Height,  orthometric  -  H  -  The  linear  distance,  measured  along 
the  local  vertical,  from  the  geoid  to  a  given  point..  Height  is 
positive  along  the  outward  normals  to  the  geopotential  surfaces. 

Height,  navigator's  -  hN  -  The  value  which  the  inertia!  navi- 
gator indicates  as  being  descriptive  of  the  third  component  of  its 
position,  in  a  direction  orthogonal  to  both  the  navigator's  lati- 
tude and  the  navigator's  longitude. 

Height  error  -  The  difference  between  the  height  and  the 
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navigator's  height.   The  two  heights  will  in  general  be  figured  at 
two  different  horizontal  locations. 
Height  error  =  hN  -  h 


CHAPTER  II 

STATIONARY  OPERATION  OF  PERFECT  INERTIAL  NAVIGATION 

DEVICES  IN  A  NON-NORMAL  FIELD 

Conceptually,  the  simplest  form  of  inertial  navigation  device 
is  one  in  which  the  reference  frame  maintains  a  fixed  orientation 
in  inertial  space,  and  which  contains  accelerometers  mounted 
along  three  mutually  perpendicular  axes.   Given  an  initial  posi- 
tion in  space,  and  an  initial  velocity,  the  accelerometer  outputs, 
when  properly  compensated  to  remove  the  effects  of  mass  attraction, 
are  doubly  integrated  to  provide,  at  any  time,  the  current  posi- 
tion.  Although  the  orientation  of  these  mutually  perpendicular 
axes  is  arbitrary,  for  ease  of  consideration  we  will  assume  that 
they  correspond  to  the  axes  of  the  Y  coordinate  system  as  de- 
scribed in  Chapter  I. 

Equations  are  developed  in  Appendix  I  to  permit  the  determi- 
nation of  acceleration  in  inertial  space  from  the  values  sensed  by 
the  errorless  accelerometers  adopted  for  this  study,  under  the 
assumption  that  gravitational  forces,  which  are  sensed  with  accel- 
erations, are  normal <»   This  is  the  assumption  used  in  currently 
operational  navigation  devices.   A  similar  development,  using 
hypothetical  accelerations  and  gravitational  field,  provides  us 
hypothetical  values  which  we  consider  to  be  sensed  by  our  error- 
less accelerometers.   After  completing  the  general  development, 
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we  consider  the  specific  case  of  a  navigator  who  is  stationary 
relative  to  the  earth  at  a  point  where  the  geoid  and  an  adopted 
reference  ellipsoid  coincide.   The  International  Ellipsoid  was 
chosen  for  convenience,  even  though  other  ellipsoids  are  presum- 
ably closer  to  the  true  size  and  shape  of  the  geoid.   Assumed 
values  of  gravity  anomaly  and  deflection  of  the  vertical  compo- 
nents were  assigned  to  the  hypothetical  point  to  which  the  navi- 
gator was  rigidly  attached,  and  these,  when  inserted  into  the 
last;  equation  of  Appendix  I,  provided  the  apparent  motion  of  the 
navigator  under  the  influence  of  non-normal  gravitation. 

Figures  2,  3»  and  k   show  the  north  geodetic  error,  the  east 
geodetic  error  and  the  height  error  caused  by  components  of  non- 
normal  gravitation.  The  term  geodetic  used  in  connection  with  the 
word  error  indicates  that  the  navigator's  position  is  equal  to  the 
geodetic  position  plus  the  geodetic  error.   Computations  were  ter- 
minated when  the  height  error  reached  +_  20  kilometers,  which  is 
well  beyond  the  limits  of  certain  assumptions  made  in  the  mathe- 
matical approximations  employed.   It  should  be  noted  that  the  ordi- 
nate scale  for  height  error  (figure  4)  is  one-tenth  that  for  the 
north  error  and  east  error,  and  that  once  the  height  error  started 
increasing,  it  was  not  a  slowly  changing  function  as  had  been 
assumed  in  obtaining  the  average  values  of  the  gravitation  func- 
tion.  The  results  do  indicate,  however,  that  this  approach  is  not 
satisfactory. 

If  we  consider  the  case  where  we  supplement  the  inertial  nav- 
igation information  with  an  independent  determination  of  the 
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17 
elevation,  we  are  able  to  limit  the  rapid  propagation  of  height 
error,  although  we  may  have  a  fairly  constant  error  of  the  uncer- 
tainty in  the  geoidal  undulation.   By  using  the  independent  deter- 
mination of  height,  we  not  only  reduce  the  height  error  to  zero, 
but  also  reduce  the  component  of  velocity  normal  to  the  ellipsoid 
to  an  independently  measured  value  at  each  computed  position.   In 
our  specific  case,  the  independently  measured  normal  velocity  is 
zero  since  the  navigator  is  fixed  to  the  earth.   The  results  of 
using  independent  height  measurements,  depicted  as  before  in  the 
form  of  north  geodetic  error  and  east  geodetic  error,  are  shown  in 
figures  5  and  6. 

From  a  casual  inspection  of  figures  5  and  6  it  would  appear 
that  there  is  a  slight  damping  of  the  north  error  where  the  deflec- 
tion of  the  vertical  is  in  the  meridian  (i.e.,  £  4   0«  tj  =0),  with 
a  time  dependent  increase  in  the  east  error,  whereas  there  is  a- 
slight  damping  with  time  of  the  east  error  when  the  deflection  of 
the  vertical  is  in  the  prime  vertical  (i.e.,  £  =  0,  ij  /  0),  with 
an  increase  with  time  in  the  north  error.   Figure  7  is  a  plot  of 
the  same  data  as  in  figures  5  and  6,  but  with  north  geodetic  error 
versus  east  geodetic  error,  showing  that  the  navigator 's  position 
is  oscillating  with  the  Schuler  period,  and  at  the  same  time  ro- 
tating a6  a  Foucault  pendulum  relative  to  the  earth.   The  example 
was  computed  for  a  latitude  of  forty  degrees,  and  the  position  of 
the  navigator  marked  at  every  fifth  minute  with  a  small  square. 
Although  the  components  of  the  navigation  error  are  time  dependent 
the  magnitude  of  the  maximum  positional  error  with  respect  to  the 
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Figure  7»   Horizontal  geodetic  error  of  indicated  position 
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astronomic  position  remains  constant.   The  magnitude  of  the  maxi- 
mum positional  error  with  respect  to  the  geodetic  position  varies 
between  twice  the  deflection  of  the  vertical  and  the  square  root 
of  two  times  the  deflection  of  the  verticalo 

The  next  step  in  conceptual  complexity  of  the  inertial  navi- 
gator, over  having  the  orientation  of  the  reference  frame  fixed  in 
inertial  space,  is  to  have  the  reference  frame  rotate  at  the  same 
rate  as  the  earth,  about  an  axis  parallel  to  the  earth's  rotation 
axis*  This  involves  only  a  constant  rotational  rate  about  one 
axis,  and  is  independent  of  the  navigator's  position  on  the  earth. 
The  output  of  the  accelerometers  is  similar  to  that  discussed 
above,  except  that  the  effect  of  the  centripetal  acceleration  and 
mass  attraction  at  a  point  maintains  a  fixed  vector  relationship 
to  the  accelerometers,  rather  than  just  maintaining  a  fixed  magni- 
tude, and  the  effects  of  Coriolis  acceleration  must  be  accounted 
for. 

Equations  are  developed  in  Appendix  II  for  the  case  of  the 
rotating  reference  frame,  analogous  to  those  developed  in  Appendix 
1  for  the  inertially  fixed  orientation  of  the  reference  frame. 
The  results  obtained  by  using  equations  II. 8  are  the  same  as  those 
obtained  from  using  equations  I. 10,  within  the  limits  of  accuracy 
of  the  mathematical  approximations  made  in  the  derivation  of  the 
equations. 

The  next  step  in  complexity  of  instrumentation  of  inertial 
navigation  systems  is  to  drive  the  platform  holding  the  acceler- 
ometers in  such  a  way  that  the  gravity  vector  always  remains 
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normal  to  the  platform.   In  such  a  system  for  surface  navigation 
we  can  eliminate  one  accelerometer ,  the  one  sensing  the  vertical 
component  of  acceleration,  and  measure  accelerations  only  in  the 
horizon  plane.   This  system  requires  that  at  every  point  where 
our  navigator  is  likely  to  be  we  know  the  detailed  shape  of  the 
geoid  so  that  we  may  drive  the  stable  platform  at  the  meridional 
local  vertical  rate  about  the  east-west  axis,  at  (earth  rate  plus 
the  astronomic  longitude  rate)  times  the  cosine  of  the  astronomic 
latitude  about  the  north-south  axis,  and  at  (earth  rate  plus  the 
astronomic  longitude  rate)  times  the  sine  of  the  astronomic  lati- 
tude about  the  vertical  axis.  This  will  provide  a  continuous 
indication  of  astronomic  latitude  and  longitude?  however,  since 
the  detailed  shape  of  the  geoid  is  known  wherever  we  navigate, 
according  to  the  requirements  of  the  system,  the  astronomic  posi- 
tion can  be  converted  to  geodetic  and  the  problem  is  solved. 

A  variety  of  things  have  conspired  to  militate  against  the 
successful  operation  of  such  a  system.   One  of  the  problems  is 
that  of  sheer  complexity  in  trying  to  describe  the  geoid,  to  com- 
pute radii  of  curvature,  and  to  drive  the  stable  platform  at  the 
correct  rates.   Another  point  which  should  not  be  overlooked  is 
our  near-total  lack  of  knowledge  of  the  shape  of  the  geoid  in  the 
ocean  areas;  however,  for  the  purpose  of  this  investigation  we 
shall  assume  that  we  do  know  the  slope  of  the  geoid  with  respect 
to  the  adopted  ellipsoid  at  our  geodetic  position. 

The  practical  implementation  of  this  type  system  has  taken 
several  different  forms,  all  of  them  assuming  that  gravity  is 
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normal  and  that  the  earth  is  an  ellipsoid,  and  that  the  stable 

platform  is  "locally  level."  That  some  of  these  assumptions  are 

not  fulfilled  is  obvious,  but  the  degree  of  lack  of  fulfillment 

may  be  illustrated  by  the  description  of  the  inertial  navigation 

equipment  used  by  the  U.S. So  Skate  in  her  transit  of  the  Arctic 

Ocean  as  given  by  Dozier: 

The  N6  platform  subsystem  mechanization  to  be  described 
may  be  regarded  as  theoretically  exact  and  error-free 
if  it  moves  on  a  perfectly  spherical  surface  and  there 
are  no  imperfections  in  the  components,  such  as  gyro- 
scopic drifts.  (2) 

Tinder  the  assumption  of  the  spherical  earth  it  was  possible 
to  divorce  platform  stabilization  from  computer  operation,  per- 
mitting the  continued  operation  of  the  system  in  the  event  of  a 
computer  malfunction,  which  is  certainly  a  very  desirable  feature; 
however,  this  system  does  introduce  certain  errors  into  the  com- 
puted position  because  of  the  fact  that  the  platform  is  not  always 
level.   Additionally,  in  other  developments,  as  a  first  order 
approximation,  height  has  traditionally  been  considered  to  lie 
along  the  direction  the  platform  is  being  maintained  perpendicular 
to,  whether  it  be  the  geocentric  radius  vector,  the  gravitation 
vector,  the  gravity  vector  or  the  normal  to  the  ellipsoid. 

In  order  to  get  away  from  the  various  approximations  made  in 
the  literature,  the  kinematic  equations  of  motion  for  an  inertial 
navigator  on  an  ellipsoid  of  revolution  are  derived  in  Appendix 
III,  being  presented  as  equations  III08.   Following  that  deriva- 
tion, a  treatment  analogous  to  that  in  Appendixes  I  and  II  for  the 
case  of  a  stationary  navigator  in  a  non-normal  field  is  carried 


out,  resulting  in  the  average  accelerations  over  a  time  step  shown 
in  equations  III.16  and  III.17<> 

When  these  equations  were  evaluated  with  the  conditions  of 
deflection  of  the  vertical  and  gravity  anomalies  used  in  the  two 
cases  above,  and  the  results  doubly  integrated,  the  errors  were 
the  same  as  before,  indicating  that  the  alignment  of  the  platform 
axes  with  the  normal  to  the  ellipsoid  at  the  navigator's  position 
rather  than  at  the  geodetic  position  was  a  source  of  error  of 
second  order,  and,  under  the  conditions  of  non-normality  of  field 
which  may  reasonably  be  expected,  may  be  disregarded o   Although 
the  value  of  position  geodetic  error  approaches  one  mile  at  the 
worst  condition,  the  difference  in  error  between  the  system  ro- 
tating with  the  earth  and  the  system  seeking  the  ellipsoidal  nor- 
mal is  a  maximum  of  0*51  meters.. 

There  are  many  variations  possible  in  the  design  of  the  sen- 
sor orientation,  but  they  are  only  minor  changes  from  the  three 
cases  considered  here*   They  may  reasonably  be  expected  to  react 
in  a  manner  similar  to  those  investigated,  so  a  detailed  investi- 
gation of  each  type  will  not  be  conducted. 


CHAPTER  III 

OPERATION  OF  A  PERFECT  INERTIAL  NAVIGATION  DEVICE 

TRANSITING  A  NON-NORMAL  FORCE  FIELD 

In  the  previous  chapter  the  operation  of  various  basic  types 
of  inertial  navigation  devices  was  investigated  for  the  restric- 
tive case  of  being  stationary  with  respect  to  the  surface  of  the 
earth.   Since  each  of  the  three  basic  types  behaved  in  the  same 
manner,  within  the  limitations  of  the  mathematical  methods  em- 
ployed, we  will  now  concentrate  on  the  type  which  is  representa- 
tive of  the  majority  of  instruments  now  available,  i.e»,  the  case 
wherein  one  axis  attempts  to  remain  aligned  with  the  normal  to 
the  reference  ellipsoid. 

Equations  are  developed  in  Appendix  IV  for  the  evaluation  of 
the  acceleration  in  latitude  and  longitude,  accepting  the  fact 
that  an  independent  means  must  be  used  to  determine  the  height  of 
the  navigator  if  this  coordinate  is  desired.   The  navigation 
problem  undertaken  using  these  equations  considered  the  case  of  a 
surface  vessel,  which  maintained  a  constant  speed  of  twenty  knots 
and  a  constant  heading  of  thirty  nine  and  one-half  degrees.   The 
fact  that  the  vessel  remained  on  the  surface  constrained  the 
orthometric  height  to  be  zero,  and  the  navigator's  estimate  of 
the  ellipsoidal  height  rate  to  be  zero. 
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A  geoidal  undulation  varying  from  zero  to  one  hundred  meters 
was  considered,  with  a  complete  undulation  cycle  occurring  in  ten 
degrees  of  latitude,  and  fourteen  and  two-tenths  degrees  of  longi- 
tude*  Although  the  height  coordinate,  measured  from  the  geoid, 
remained  zero,  the  height  rate  with  respect  to  the  ellipsoid  was 
not  zero,  because  of  the  geoidal  undulationso   Figures  8  through 
22  show  the  values  of  the  components  of  the  deflection  of  the 
vertical,  £,  and  t|,  as  a  function  of  time.  Figures  8,  9,  and  10 
also  show  the  error  of  position,  relative  to  the  astronomic  posi- 
tion.  In  this  simulation,  which  ran  over  a  navigation  period  of 
39  hours,  the  maximum  error  relative  to  the  astronomic  position 
was  5*6  meters,  occurring  just  after  the  twenty-fifth  hour.   This 
figure  can  be  compared  to  the  error  with  respect  to  the  geodetic 
position  of  as  much  as  1750  meters  (0.95  miles). 

The  results  of  this  test  can  be  generalized  as  indicating 
that,  while  assuming  that  navigation  is  taking  place  on  an  ellips- 
oid and  that  the  force  of  gravity  is  everywhere  normal  to  this 
surface,  the  inertial  navigator  (mechanically  perfect),  when  ini- 
tially set  to  proper  position  and  velocity,  will  track  the  actual 
direction  of  gravity  on  the  geoidal  surface  to  a  very  high  degree 
of  accuracy,  and  will  provide  the  angle  between  the  gravity  vector 
and  the  equatorial  plane  as  latitude  and  the  angle  between  the 
plane  which  contains  the  gravity  vector  and  which  is  parallel  to 
the  rotational  axis  of  the  earth  and  the  Greenwich  meridional 
plane  as  the  longitude.   Thus,  if  a  map  of  the  deflection  compo- 
nents were  available  with  an  argument  of  astronomic  position,  it 
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would  be  a  simple  matter  to  transform  the  position  so  obtained  to 
the  geodetic  position.   The  effect  of  changing  the  gravity  anomaly 
from  plus  to  minus  was  investigated  for  a  period  of  1000  minutes, 
with  a  maximum  difference  of  position  of  eleven  centimeters  which 
occurred  at  825  minutes  when  the  value  of  the  anomaly  was  +  46.3 
milligals.   Although  the  anomaly  value  went  as  high  as  +_  71° 6 
milligals  during  this  simulation,  the  position  differences  between 
the  two  conditions  decreased  after  825  minutes. 

Although  this  study  considers  only  a  perfect  inertial  system 
and  blithely  ignores  all  the  problems  that  have  been  plaguing 
engineers  for  years,  we  know  that  position  errors  do  creep  into 
the  inertial  navigator  over  a  period  of  time.   For  this  reason,  it 
is  necessary  to  reset  the  equipment  periodically  to  the  correct 
position  and  velocity.   Let  us  assume  that  some  external  naviga- 
tion system,  such  as  the  Navy  Navigational  Satellite  System  de- 
scribed by  Russell  (9),  will  provide  errorless  positioning  infor- 
mation in  an  absolute  system,  and  additionally  will  provide  error- 
less information  on  the  component  of  velocity  in  the  north  direc- 
tion. We  will  use  this  external  source  of  errorless  information 
every  eight  hours  for  resetting  our  inertial  navigator  to  geodetic 
position  and  actual  velocity  north.   The  results,  figures  11,  12 
and  13,  indicate  that  a  positional  error  of  3*2  kilometers  (1.73 
miles)  has  been  reached,  with  respect  to  the  geodetic  position, 
within  A-5  minutes  of  the  second  reset.   In  view  of  the  fact  that 
the  navigator  attempts  to  track  the  local  vertical,  let  us  consid- 
er the  error  with  respect  to  the  astronomic  position,  but  under 
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the  same  resetting  conditions.   As  shown  in  figures  l*f,  15  and  16, 
the  maximum  error  is  reduced  to  1.66  kilometers  (0.89  miles), 
which  is  still  considered  excessive  for  any  type  of  survey,  or  for 
any  other  purpose  requiring  precise  navigation. 

In  order  to  exploit  the  fact  that  the  inertial  navigator 
seeks  the  local  vertical,  a  simulated  run  was  made  wherein  the 
position  was  reset  not  to  the  geodetic  position,  but  to  the  astro- 
nomic position,  while  the  velocity  north,  as  obtained  without  er- 
ror from  the  Navy  Navigational  Satellite  System,  was  used  without 
change.   The  results  are  shown  in  figures  17,  18  and  19.   The 
reader  is  cautioned  about  the  change  of  ordinate  scale  between  the 
different  figures.   Although  figures  17,  18  and  19  look,  at  first 
glance,  as  discouraging  as  figures  11  through  16,  the  maximum 
error,  with  respect  to  the  astronomic  position,  is  less  than  52 
meters  (0.028  miles).   This  error  is  of  a  reasonable  magnitude,  and 
should  be  acceptable  for  most  purposes;  however,  for  purposes  of 
special  surveys  there  are  instances  in  which  this  error  may 
be  excessive. 

In  navigating  through  an  area  where  the  deflection  of  the 
vertical  is  changing,  the  astronomic  latitude  rate  of  the  naviga- 
tor consists  of  two  independent  components.   One  of  these  is,  of 
course,  the  velocity  of  the  vessel  over  the  surface  of  the  ellips- 
oid, which,  when  divided  by  the  appropriate  radius  of  curvature, 
becomes  the  angular  rate  of  the  vessel.   This  item,  for  the  lati- 
tude rate,  was  accounted  for  in  the  analysis  above.   The  second 
item  which  must  be  considered  is  the  time  rate  of  change  of  the 
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deflection  components.,   In  the  present  case,  in  which  we  consider 
external  data  affecting  only  position  and  latitude  rate,  we  set 

*     n   •     •     • 

where 

<p4   is  the  astronomic  latitude  rate  of  the  navigator 

Vn  is  the  component  of  velocity  in  the  north  direction  as 
provided  by  an  external  data  source 

M,j   is  the  meridional  radius  of  curvature  for  the  navigator's 
latitude 

4,   is  the  time  rate  of  change  of  the  meridional  component 
of  the  deflection  of  the  vertical. 

Figures  20,  21  and  22  show  the  results  of  this  simulation. 
The  maximum  error  introduced  in  this  case,  with  respect  to  the 
astronomic  position,  was  less  than  12  meters  (O.OO65  miles),  which 
is  considered  completely  acceptable  for  even  the  most  stringent 
navigation  requirements.   It  should  be  noted  that  the  geoidal  un- 
dulation was  not  considered  in  determining  the  astronomic  latitude 
rate  of  the  navigator.   Presumably  better  results  could  have  been 
obtained  if  this  had  been  accounted  for,  but  since  the  magnitude  of 
the  undulation  is  so  small  compared  to  the  ellipsoidal  radius  of 
curvature,  which  ranges  between  6.3  and  6.4  megameters,  it  was 
considered  an  unnecessary  complication  to  include  it  here. 

The  environment  in  this  test  case  is  considered  to  be  quite 
severe,  in  that  an  undulation  of  100  meters  is  periodic  in  less 
than  850  miles.   For  a  specific  case  of  navigation,  considering 
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the  slope  of  the  geoid  in  the  area  of  interest  and  the  maximum 
speed  of  the  vessel  to  be  navigated,  it  may  well  be  that  the  in- 
elusion  of  £  is  not  warranted  and  that  the  velocity  north  divided 
by  the  meridional  radius  of  curvature  is  a  sufficiently  accurate, 
approximation  of  the  astronomic  latitude  rate  to  provide  the  navi- 
gational accuracy  requiredo   In  any  event,  it  is  considered  that 
if  5  is  to  be  considered,  it  should  be  determined  for  the  geodetic 
position  from  a  plot  of  deflection  components,  and  entered  man- 
ually into  the  navigational  computer  rather  than  utilizing  storage 
in  the  computer  to  describe  the  geoid  and  to  compute  the  deflec- 
tion rate  within  the  computer* 

We  have  considered  only  the  use  of  the  Navy  Navigational 
Satellite  System  to  provide  external  data  for  the  inertial  navi- 
gator.  Other  systems,  such  as  bottom  navigation,  could  provide 
not  only  the  information  considered  above,  but  also  the  velocity 
in  the  east  direction,  which  would  be  handled  in  a  manner  analo- 
gous to  the  velocity  north  to  further  reduce  the  error  after  reset 
of  the  inertial  navigator . 

In  order  to  navigate  and  obtain  the  geodetic  position  at  all 
times,  it  would  be  necessary  to  compute  the  components  of  gravity 
in  the  horizontal  plane  at  each  step  of  the  navigational  position 
updating  which  takes  place  several  times  a  second  in  most  naviga- 
tors.  This  would  require  extensive  computer  time  and  increase  the 
storage  requirements  greatly.   By  the  simple  expedient  of  reset- 
ting the  navigator's  position  to  the  astronomic  rather  than  to  the 
geodetic  position,  and  by  inserting  the  astronomic  latitude  and 
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longitude  rates  (which  include  the  rates  of  change  of  the  deflec- 
tion components),  we  are  able  to  accomplish  precise  navigation 
without  introducing  an  increase  in  the  time  required  per  itera- 
tion, or  increasing  the  size  of  the  navigational  computer..   In 
order  to  use  the  information  thus  produced  by  the  navigator  it  is 
necessary  to  consider  the  difference  between  the  local  vertical 
and  the  ellipsoidal  normal,  but  this  can  be  done  by  graphical 
means  or  in  another  computer,  and  done  only  when  required,  rather 
than  at  each  step  of  the  iterative  navigational  problem. 

We  have  considered  only  the  case  where  the  externally  sup- 
plied information  is  errorless.   Considering  the  pendulum  type 
effect  of  positional  errors,  we  can  conclude  that  an  error  in  the 
position  to  which  the  inertial  navigator  is  reset,  regardless  of 
cause,  will  cause  the  half -amplitude  of  the  error  shown  in  the 
figures  (17  through  22)  to  be  increased  by  the  amount  of  the 
error.   Since  the  Foucault  rotation  must  be  considered,  the  half- 
amplitude  of  the  individual  components  cannot  be  increased  accord- 
ing to  the  north  error  or  east  error,  as  appropriate,  but  the  mag- 
nitude of  the  horizontal  error  must  be  added  to  the  magnitude  of 
the  indicated  error.   The  positional  error  at  reset  can  be  caused 
by  either  incorrect  data  from  the  external  navigation  device,  or 
by  errors  in  the  determination  of  the  deflection  of  the  vertical. 
The  error  which  is  indicated  in  figures  17  through  22  is  the 
result  of  velocity  errors  at  the  time  of  reset,  primarily,  and  the 
result  of  the  difference  between  the  mathematical  model  we  have 
used  for  the  earth  and  the  true  earth  shape  secondarily.   If  we 
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consider  the  maximum  error  shown  in  figure  21  (12  meters)  to  be 
caused  entirely  by  a  velocity  error  at  reset,  it  represents  0.03 
knots  error.   An  error  of  0.1  knots  (0.0515  meters  per  second) 
would  result  in  a  positional  error  of  41.5  meters  half-amplitude, 
with  a  Schuler  period. 


CHAPTER  IV 

MAGNITUDE  OF  DEFLECTION  VALUES  WHICH  MIGHT 

REASONABLY  BE  ENCOUNTERED  AT  SEA 

Now  that  we  have  seen  the  effect  of  the  deflection  of  the 
vertical  on  the  operation  of  a  perfect  inertial  navigator,  we  may 
logically  turn  to  the  question  of  the  magnitude  of  deflection  of 
the  vertical  values  which  might  reasonably  be  encountered  at  sea*. 
If  a  complete  gravity  survey  of  the  world  were  available  it  would! 
be  possible  to  map  the  deflections  for  any  area  of  interest  to  see 
whether  or  not  they  were  large  enough  to  be  concerned  about,  and 
to  make  compensation  for  when  resetting  the  inertial  navigator* 
In  the  absence  of  an  adequate  world-wide  survey,  it  is  necessary 
to  turn  to  other  means  to  determine  whether  or  not  we  need  an  ac- 
curate determination  of  the  deflection  of  the  vertical  in  the 
ocean  areas* 

Rice  (8),  in  his  study  of  the  area  300  kilometers  by  1800 
kilometers  reaching  from  Cape  Kennedy,  Florida,  to  the  southwest 
corner  of  Puerto  Rico,  which  is  largely  an  ocean  area,  obtained 
deflections  as  great  as  56.7"*   This  area  can  be  considered  as 
extreme,  since  on  one  side  of  it  there  are  extensive  banks  which 
are  only  a  few  fathoms  deep,  with  the  other  side  plunging  rapidly 
to  several  thousand  fathoms,  and  we  would  be  inclined  to  avoid  it 
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if  we  wanted  to  use  inertial  navigation  without  compensating  for 
the  deflection  of  the  vertical.   Needless  to  say,  such  a  value  is 
not  typical  of  the  entire  area. 

Let  us  consider  bathymetric  features  which  are  covered  by  a 
minimum  of  1000  fathoms  of  water,  so  that  the  first  1000  fathoms 
below  the  navigating  unit  is  homogeneous*  There  are  features  in 
the  ocean  areas  which  are  of  such  an  extent  that  they  may  be  con- 
sidered to  be  of  infinite  length,  and  may  be  approximated  as  being 
made  of  a  combination  of  several  bodies  which  are  rectangular  in 
cross  section.   Heiland  (3)  has  given  formulas  to  compute  the 
attraction  of  infinite  rectangles  as  follows: 
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where  a   is  the  density  difference  between  the  material  of  the  dis- 
turbing mass  and  that  of  the  homogeneous  material  surrounding  it. 


50 

The  above  equation  may  be  rewritten  for  the  horizontal  rather 
than  the  vertical  component  of  the  change  in  gravity  created  by 
the  disturbing  mass,  in  that  horizontal  direction  perpendicular  to 
the  axis  of  the  disturbing  mass,  as 
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We  will  assume  that  we  are  working  on  a  plane,  so  that  Agx  is 

— ♦ 
perpendicular  to  g,  making  it  possible  to  represent  the  deflection 

by  — *"-.  Three  cases  have  been  investigated  which  are  considered* 

O 

as  typical  of  the  types  of  bathymetric  features  which  may  be  found 
in  mid  ocean,  and  near  which  operations  may  presumably  be  con- 
ducted. 

CASE  I.   The  Mid-Atlantic  Ridge  has  been  approximated  by 
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Figure  23 «>   Approximation  of  Mid-Atlantic  Ridge 

three  blocks  of  500  fathom  thickness,  with  widths  of  70,  100  and 
130  miles,  spanning  the  depths  between  1000  and  2500  fathoms,  as 
shown  in  figure  23 • 
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CASE  II.   The  ridge  southeast  of  the  Japanese  Islands,  near 
<p  =  32  ,  \.   =  139  1  is  similar  to  the  Mid-Atlantic  Ridge,  but  rises 
out  of  a  more  shallow  ocean.   It  has  been  approximated  by  two 
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Figure  2k.      Approximation  of  ridge  south  of  Japan 

blocks  of  500  fathom  thickness,  with  widths  of  77«5  and  112.5 
miles,  spanning  the  depths  between  1000  and  2000  fathoms,  as 
shown  in  figure  2k. 

CASE  III.   The  northern  part  of  the  Marianas  Trench,  near 
<p  =  32°,  \.   =  142°,  drops  out  of  an  ocean  3500  fathoms  deep  to  a 
depth  of  more  than  5000  fathoms.   It  has  been  approximated  by 
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Figure  25.   Approximation  of  Marianas  Trench 

three  blocks,  each  500  fathoms  thick,  with  widths  of  43,  27  and 
13  miles,  spanning  the  depths  from  3500  to  5000  fathoms,  as  shown 
in  figure  25.   It  is  by  far  the  smallest  feature  of  the  three 
considered,  and  is  covered  by  three  and  one-half  times  as  thick  a 
layer  of  homogeneous  matter;  therefore,  we  might  logically  expect 


52 

a  much  smaller  deflection  in  Case  III  than  in  Cases  I  and  II. 

The  magnitude  of  the  deflection  in  each  case  is  shown  in 
figures  26,  27  and  28,  as  a  function  of  the  distance  from  the 
centerline  of  the  disturbing  masses..   The  maximum  deflection  is, 
of  course,  from  Case  I,  being  32.9"  at  50  miles  from  the  center- 
line,  with  a  deflection  of  over  11"  remaining  2  therefrom.   In 
Case  II  a  maximum  deflection  of  2k.7u   was  found  56  miles  from  the 
centerline,  with  over  7"  remaining  2  from  the  disturbance.   Case 
III,  with  the  minimum  effect,  had  a  maximum  deflection  of  18.7" 
at  20  miles  from  the  centerline,  diminishing  to  2.9"  at  a  distance 
of  2°. 

Ordinarily,  a  deflection  value  of  less  that  3"  of  arc  would 
be  considered  negligible;  however,  if  we  reset  the  inertial  navi- 
gator to  the  true  position  without  considering  the  deflection  of 
the  vertical  (i.e.,  to  the  geodetic  position),  after  kZ   minutes 
we  would  have  a  position  geodetic  error  of  0.1  miles.   The  im- 
portance of  this  error  would  depend  on  the  mission  of  the  vessel 
using  the  navigator,  and,  while  perfectly  acceptable  for  many 
missions,  this  error,  in  conjunction  with  the  error  build-up  in 
practical  navigators,  could  be  the  cause  of  the  failure  of  some 
other  mission. 

It  should  be  noted  that  the  bathymetry  used  in  Case  II  and 
Case  III  occur  with  a  longitude  difference  of  only  3  ,  or  a  spher- 
ical separation  of  about  2.5° •   If  the  navigator  were  between  the 
two  disturbances,  the  effects  would  be  additive  in  magnitude.   In 
each  of  these  cases  we  have  considered  the  effect  of  only  one 
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disturbing  feature,  and  have  completely  ignored  the  compensating 
effect  of  any  isostatic  compensation,  which  would  be  of  increasing 
importance  as  the  distance  from  the  disturbance  increased. 

If  we  consider  the  world  as  a  whole,  rather  than  isolated 
features,  we  may  work  with  the  gravitational  potential  as  deter- 
mined from  satellites,  and  see  what  the  undulations  of  an  equi- 
potential  surface  derived  from  satellite  observations  are  with 
respect  to  an  ellipsoid  of  rotation.   First,  let  us  recognize  that 
the  potential  determined  from  satellites  at  altitudes  of  one  or 
more  megameters  cannot  adequately  describe  the  potential  on  the 
surface  of  the  earth,  so  the  undulations  determined  are  not  those 
of  the  geoidc   They  represent,  rather,  the  undulations  of  an  equi- 
potential  surface  of  a  harmonic  potential  function  in  which  the 
coefficients  of  the  higher  order  terms  have  been  inadequately 
determined  because  of  the  attenuating  effect  of  the  distance  be- 
tween the  masses  creating  the  potential  and  the  altitude  at  which 
the  effects  are  measured,  in  conjunction  with  the  inadequacy  of 
the  tracking  equipment  to  detect  small  movements.   Additionally, 
the  insufficiency  of  numbers  of  tracking  equipments  to  cover  the 
entire  trajectories  of  these  satellites  gives  a  fictitious  poten- 
tial function. 

In  spite  of  the  shortcomings  of  the  potential  function  which 
we  obtain  from  satellite  observations  when  we  want  to  apply  that 
function  at  the  surface  of  the  earth,  the  undulations  obtained 
from  the  satellite  derived  potential  should  correspond  in  a  gross 
way  to  those  obtained  from  the  gravimetric  method,  and  the  slope 
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of  the  undulating  surface  with  respect  to  an  approximating  ellips- 
oid should  agree  with  the  average  deflections  of  the  vertical  ob- 
tained by  the  gravimetric  method  when  compared  in  an  area  of  great 
enough  extent,  as  determined  by  the  degree  to  which  the  satellite 
potential  is  accurately  computed..   Thus,  if  we  were  to  compute  the 
slope  of  the  satellite  potential  surface  relative  to  the  ellipsoid, 
we  would  not  expect  to  agree  in  any  measure  with  the  deflection  of 
the  vertical  as  obtained  from  gravimetry?  however,  if  we  obtain 
the  average  slope  of  this  surface  over  an  extended  distance,  we 
would  expect  fair  agreement  with  the  average  of  the  deflections  of 
the  vertical  over  the  same  distance  <> 

The  gravitational  field  representation  developed  by  Anderle 
(1),  which  is  developed  through  the  seventh  degree,  sixth  order, 
was  converted  to  undulations  as  outlined  by  Nicolaides  and  Macom- 
ber  (7).   The  formula  for  the  undulations  is 


r  7    n 


C  »  P 


£    £  P/Csin*: 


C      cosmX   +  S      sinm\. 
nm  nm 


n=2     m»0 

-   J3Pa  (sintp)    -  J4P4  (simp) 

where 

p  is  the  radius  vector 

P   are  Legendre  Polynomials 

P  m  are  Associate  Legendre  Functions 
n 

C   and  S   are  coefficients  in  the  gravitational  potential 
nm      nm 

expansion,  and 
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Ja  and  J4  apply  to  the  reference  ellipsoid 

The  Js  and  J4  used  were  for  a  reference  ellipsoid  flattening 
of  1/298.3°   From  these  undulations,  deflections  were  computed 
from  the  equations 


1  = 


M  A<p 
AC 


N  coscp  A\ 


The  results  of  this  analysis  showed  that  in  the  North  Atlantic 
Ocean,  in  a  region  bounded  on  the  south  by  <p  =  30  ,  between  the 
meridians  of  -  30  and  -  kQ    ,  and  on  the  north  by  cp  =  70°,  between 
the  meridians  of  -  30  and  -  80  ,  the  average  slope  is  on  the  order 
of  seven  and  one-half  seconds  of  arc.   This  is,  admittedly,  the 
largest  ocean  area  which  shows  a  systematic  slope  of  this  order  of 
magnitude.   If  a  perfect  inertial  navigator  were  reset  to  the  geo^- 
detic  position  at  a  spot  where  the  average  deflection  were  appli- 
cable, the  position  geodetic  error  after  a  period  of  kZ   minutes 
would  be  one-quarter  mile.   Since  the  average  error  is  of  the 
order  of  seven  and  one-half  seconds,  it  is  reasonable  to  expect 
some  areas  where  the  deflection  is  larger,  resulting  in  position 
geodetic  errors  larger  than  one-fourth  mile  if  the  navigator  were 
reset  without  compensating  for  the  effect  of  the  deflection  of  the 
local  vertical. 

There  is  nothing  in  this  chapter  which  should  be  interpreted 
to  indicate  that  deflections  of  the  vertical  of  any  given  magni- 
tude do  exist  in  the  ocean  areas.   The  purpose  of  this  chapter  is 
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to  show  that  one  cannot  reasonably  exclude  the  possibility  that 
deflections  exist  in  deep  ocean  areas,  which  are  of  such  a  magni- 
tude that  they  could,  if  ignored,  introduce  sufficient  error  into 
an  inertial  navigator  to  cause  the  failure  of  some  mission  which 
requires  precise  navigation  for  its  successful  accomplishment o 


CHAPTER  V 
CONCLUSIONS  AND  RECOMMENDATIONS 

A  mechanically  perfect  inertial  navigator  is  extremely  sensi- 
tive to  the  effects  of  the  magnitude  of  the  gravity  vector  if  an 
attempt  is  made  to  navigate  along  that  vector .   A  slight  error  in 
the  magnitude  of  gravity  will  cause  an  unbounded  error  to  develop 
in  the  vertical  direction..   This  unbounded  error  source  makes  the 
inertial  navigator  totally  unsuited  for  the  navigation  task  in 
this  direction.   If  we  consider  a  case  in  which  there  is  an  in- 
dependent input  for  height  information,  and  the  inertial  navigator 
is  used  only  for  navigation  in  the  horizontal  plane,  the  magnitude 
of  gravity  is  of  no  practical  importance. 

When  we  use  a  mathematical  model  for  our  navigation  problem 
which  is  based  on  an  ellipsoid  of  rotation  for  the  figure  of  the 
earth,  with  position  defined  by  the  normal  to  that  ellipsoid,  we 
find  that  the  inertial  navigator  seeks  the  local  vertical  rather 
than  the  normal  to  the  fictitious  reference  ellipsoid.   If  we  set 
the  navigator  to  the  geodetic  position,  the  navigator's  position 
will  oscillate  through  twice  the  angle  of  the  deflection  of  the 
vertical,  with  a  period  of  a  Schuler  pendulum,  and  the  navigator's 
position  will  lie  near  a  plane  rotating  about  the  astronomic 
position  with  the  angular  velocity  of  a  Foucault  pendulum  relative 

to  the  earth. 
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If  we  had  a  perfect  inertial  navigator,  and  set  the  position 
initially  to  the  astronomic  position  while  stationary  with  respect 
to  the  earth,  then  the  navigator  would  provide  us  continually  with 
positions  within  a  few  meters  of  the  astronomic  positions  as  we 
traversed;  an  area  with  varying  deflections  of  the  vertical- 

Unfortunately,  a  mechanically  perfect  inertial  navigator  has 
not  been  developed,  so  that  errors  do  creep  into  the  system,  mak- 
ing it  necessary  to  reset  the  navigator  periodically.   Resetting 
to  the  geodetic  position  is  unsatisfactory  since  the  navigator's 
position  will  then  oscillate  with  twice  the  amplitude  of  the  de- 
flection of  the  vertical.   In  order  to  avoid  this  oscillation  and 
at  the  same  time  retain  the  current  level  of  sophistication  in  the 
computer  programs,  it  is  mandatory  that  the  position  for  resetting 
is  that  defined  as  the  astronomic  position,  and  that  the  angular 
velocities  we  accept  at  the  time  of  resetting  be  the  astronomic 
latitude  rate  and  the  astronomic  longitude  rate. 

If  accurate  navigation  is  required  through  the  use  of  an 
inertial  navigator,  it  is  necessary  that  some  means  be  available 
for  resetting  position  to  that  of  the  astronomic  position  and 
angular  velocity  to  the  astronomic  rates,  or,  alternatively,  to 
know  the  gravity  vector  precisely,  and  to  compensate  for  the 
geodetically  horizontal  components  if  the  navigator  is  reset  to 
the  geodetic  position  and  the  angular  velocity  is  reset  to  the 
geodetic  latitude  and  longitude  rates.   The  alternative  would 
necessitate  a  great  increase  in  the  computation  load  for  deter- 
mining positions. 
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The  ideal  case  for  resetting  the  navigator,  from  the  stand- 
point of  simplicity,  is  to  do  so  in  an  area  where  both  the  deflec- 
tion of  the  vertical  and  the  time  rate  of  change  of  the  local 
vertical  are  zero*   It  is  realized  that  this  condition  may  be 
hard  to  achieve  operationally., 

It  is  recommended  that,  where  a  valid  requirement  exists  for 
precise  navigation,  and  where  an  inert ial  navigator  is  to  satisfy 
that  requirement,  that  a  complete  gravimetric  survey  of  the  area 
be  made,  and  that  charts  showing  the  deflection  of  the  vertical  be 
prepared,  with  entering  arguments  of  geodetic  position  in  one  case, 
and  of  astronomic  position  in  the  other.   This  will  permit  taking 
a  geodetic  position  and  converting  it  to  astronomic  when  resetting 
the  navigator,  and  to  compute  the  time  rate  of  change  of  the  com- 
ponents of  the  deflection  of  the  vertical  which  are  necessary  in 
determining  the  astronomic  latitude  and  longitude  rates,  and  then 
taking  an  updated  navigator's  position  and  converting  it  to  an 
approximate  geodetic  position  when  needed. 

The  fact  that  the  effects  of  non-normal  gravity  were  consid- 
ered on  a  mechanically  perfect  inertial  navigator  should  not  be 
construed  to  indicate  that  the  writer  feels  that  the  mechanical 
imperfections  are  unimportant  and  can  be  ignored.   Rather,  he 
feels  incompetent  to  consider  them  without  extensive  collabora- 
tion with  engineers  specializing  in  the  field.   Now  that  the 
effects  of  non-normal  gravity  have  been  investigated  while  being 
isolated,  it  would  be  most  appropriate  to  undertake  a  study  com- 
bining the  effects  of  non-normal  gravity  and  mechanical  imper- 
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fections,  so  that  it  may  be  seen  where  the  influence  of  one  over- 
shadows the  influence  of  the  other. 

Although  the  underlying  assumption  throughout  this  study  has 
been  that  information  on  the  deflection  of  the  vertical  would  be 
obtained  through  gravimetric  studies,  one  should  not  rule  out  the 
possibility  of  developing  a  sensor  which  would  be  able  to  distin- 
guish between  the  mass  attraction  effects,  the  centrifugal  force, 
and  the  reaction  to  acceleration.   Such  a  sensor  would  change 
these  conclusions  and  recommendations  completely.   Similarly,  if 
a  sensor  could  be  developed  which  would  continuously  analyze  the 
gradients  of  gravity  to  carry  either  absolute  or  relative  deflec- 
tions with  the  navigating  vessel,  the  recommendations  would  have 
to  be  modified  accordingly. 


APPENDIX  I 

DEVELOPMENT  OF  FORMULAS  FOR  SENSORS 

FIXED  IN  INERTIAL  SPACE 

The  Y  coordinate  system  will  be  used  in  this  case.   The 
accelerometers,  which  sense  along  the  coordinate  axes,  respond 
to  mass  attraction  and  accelerations-.   Streeter  (10)  has  shown 
that  accelerations  due  to  rotation  about  the  sun  can  be  neglected 
and  that  the  mass  attractions  of  the  matter  outside  the  earth  are 
negligible.   Letting  As  represent  the  value  sensed  by  the  accel- 
erometers, we  may  write 

A8  *  A*,  -  G  1,1 

=  AmB  +  Ag  ,  -  G 

where  AN,  =  the  acceleration  of  the  navigator  with  respect  to 
inertial  space, 
ANB  =  the  acceleration  of  the  navigator  with  respect  to 

point  B,  said  point  being  fixed  to  the  earth  on  the 
surface  of  the  reference  ellipsoid,  and 
Ae  ,  =  the  acceleration  of  point  B  with  respect  to  inertial 
space. 
Since  point  B  is  fixed  to  the  reference  ellipsoid,  the 
acceleration  is  merely  centripetal? 
Ae,  =  Ox(QxR)  =  nxCflxNe  ) 

6^ 
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where  Nb  is  the  vector  from  the  spin  axis  to  point  B which  is  nor- 
mal to  the  reference  ellipsoid  at  point  B. 

As  =  Ox(QxN8  )  +  ANB  -  G  1.2 

In  the  problem  of  practical  navigation,  the  only  quantity- 
measured  is  the  sensed  value  provided  by  the  accelerometer . 
Since  the  value  for  the  gravitational  vector  is  imperfectly 
known,  "ft  is  approximated  by  "normal"  gravitation,  providing, 
in  general,  a  false  indication  of  the  acceleration  of  the  navi- 
gator, R*B  ° 

As  =  QxCQxHe  )  +  RnB  -  % 

Rnb  =  Anb  -  G  +  GN  1.3 

Normal  gravitation  has  been  developed  (Hirvonen  and  Moritz, 
(5))  in  terms  of  the  geocentric  radius,  R,  the  geocentric  lati- 
tude, <p,  and  the  longitude,  \,  once  the  defining  parameters  of 
the  gravitational  reference  ellipsoid  are  known.   For  the  Inter- 
national Ellipsoid  with  the  International  normal  formula  of 
gravity,  when  R  is  expressed  in  kilometers  and  G  in  gals,  we 
obtain,  the  following: 

G_  =  -  1000  F  +  0.835  888  F3  -  2.507  664  5*0032^  + 
+  0.005  118  F3  +  0.005  937  F*cos2Jp  - 

-  0.013  778  F3cosa2Jp  - 

-  0.000  007  F4    +  0.000  051*  F*cos2<p  + 

+   0.000   023  F*cos32(p   -   0.000  085  F4  cos3  2cp, 

where  F   =   3986    >290q^,      and   a     is   the    component    of   normal   gravi- 
da K 
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tation  along  the  geocentric  radius  vector.      The   component   in  the 
meridional  plane   orthogonal   to   the   radius  vector,    G- .    is 

G-   =   sinq>   cosip      -   3.3^3  351  F3    + 

+   0.003  1^9  F3    -   0.022  043   F3cos2^   + 

+   0.000  031  F*    +   0.000  026  F*cos2^   - 

-   0.000  1^5  F4cosa2^ 

The  components  of  gravitation  along  the  earth-fixed  coordi- 
nate axes  may  then  be  obtained  from 

QyO    =      GR  coscp  -   G-   sin<p|   cos\ 

G-jj-i    =    [G_   coscp   -   G-   sinip     sin\. 

Gy3    =      Gg  sinip  +   G-   cosip 

The  value  of  the  gravitation  vector  actually  existing  at  a 
given  point  is  the  sum  of  the  normal  gravitation  vector  and  the 
gravity  disturbance  vector,  6,  at  that  point.  The  gravity  dis- 
turbance vector  is  computed  as  shown  in  Hirvonen  and  Moritz  (5)» 
and  transformed  to  the  coordinate  axes  in  the  same  way  as  the 
normal  gravitation  vector. 

After  we  have  added  the  components  of  the  normal  gravitation, 
G,  and  the  components  of  the  gravity  disturbance  vector,  6 ,  it  is 
necessary  to  transform  them  from  the  earth-fixed  coordinate  system, 
X,  to  the  inertial  coordinate  system,  Y,  by  considering  the  angle 
through  which  the  earth  has  turned,  CH»  We  then  put  the  trans- 
formed components  into  equation  1.3  to  obtain  the  value  which  the 
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inertial  navigator  interprets  to  be  his  acceleration  in  inertial 
space,  RMB  . 


inb    = 


"i°    +    (Qyo    +  6^    -   Gxo)cosnt-    (Gxl    +  6xl    -   Gxl  )sinntli 


[y1    +    (G^o    +  Sjjo    -   Gxo)sinf2t+    (Q^    +  6xi    -  Gxi)cosfit]3   + 


>lmk 


[y3   ♦  ( 


Gx3    +6 


X3 


-v] 


In  the  present  case,  wherein  the  vehicle  to  be  navigated  is  a 
surface  vessel  and  remains  close  to  a  known  equipotential  surface 
rather  than  navigating  freely  in  three  dimensions,  we  may  simplify 
the  procedure  for  computing  both  the  normal  and  the  actual  gravi- 
tation by  using  gravity  anomalies,  Ag,  rather  than  gravity  dis- 
turbances, 6go   The  magnitude  of  the  actual  gravity  on  the  geo- 
potential  surface  W  =  Wo  is  normal  gravity,  computed  on  the 
spheropotential  surface  U  =  W0 ,  plus  the  gravity  anomaly.   The 
direction  in  which  it  acts  is  defined  by  the  angles  <p  +  4  and 

\.    +     Tf/C-OSipo 


Then 


G  =  g  +  C)x(ftxR) 

"&N    =  -(«yM    +  ft3  (NN    +  hN))cos(pN    cos(,Yn    +  fit)      i 

-(Yn    +Qa(NN    +   hN))costpM    sin(\N    +  fit)      j 

-  Vn    sin<pN   k 
Similarly,    we  may  write   the   expression   for   the  actual  gravi- 
tation and   account    for    the   non-normality   of   fieldo 
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G   =   -    (v   +  Ag)cos(q>   +    £)    cos(V   +    n/cos<p   +  Qt)    + 


fJa  (N  +    h)cos<p   cos(\.   +  Clt) 


(y   +  Ag)cos(<p   +   £.)    sin(\.  +   Tj/cosip  +  Clt)    + 


fla  (N  +   h)costp   sin(\   +  Clt) 


-    (V   +  Ag)sin((p  +   6.)     k 


Substituting  these  expressions  for  the  actual  and  the  normal 
gravitation  into  equation  1.3 »  w®  obtain  an  expression  for  the 
value  of  the  acceleration  used  in  the  navigation  problem. 


Rnb    = 


1°    +    (y  +  Ag)cos(cp  +   4)    cosCV  +   ri/coscp  +  Clt)   + 


+  fla  (N  +  h)coscp  cosCV  +  Qt)   - 


-   (\N    +  CI*  (NN    +  hH))cos(pN    cos(\*    +  Clt) 


i   + 


+IT1    +    (y  +  Ag)cos((p  +   O   sin(k  +  Ti/cos«p  +  Clt)   + 


>1.5 


3   + 


+  na(N  +  h)cos<p  sinCX.  +  fit)   - 
-    (yN    +  na  (NN    +  hN))co8<pN    sin(\«    +  CH) 
Y8    +    (v   +  Ag)sin(<p  +   O    -  YMsintpH 

The  gravity  formula  used  in  this  study  is  the  formula  for 
the  International  Ellipsoid,  modified  to  include  a  term  accounting 
for  free-air  reduction  where  the  height  is  positive,  and  a  term 
accounting  for  the  free-air  reduction  plus  the  effect  of  salt 
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water  in  those  cases  where  the  height  is  negative.   The  formulas 
used  are: 

a.  For  h  positive  - 

Y  *  978.049   (1  +  0.0052884  sin3<p  -  0.0000059  sina  2(p)   - 

-  0.3086  h  x  10"3  centimeters/sec3 

b.  Dor  h  negative  - 

y  =  978o 049   (1  +  0.0052884  sin3^  -  0.0000059  sina2(p)   - 

-  0.2223  h  x  10-3    centimeters/sec3 

The  determination  of  position  in  6uch  a  rectangular  system 
is  very  straight-forward,  and  may  be  developed  in  the  following 
manner. 

S=^)+^,t+     /         /      Rdudv 
0        0 

3   =  %    +  RoNB  t    +   RoB  ,  t   + 

1 

t  V  t  V 

+      /         /      RMB    du  dv  +      /         /       Rg  ,     du  dv 
0        0  0        0 

Introducing  the  components  of  the  initial  value  of  the  radius 
vector,  the  initial  velocity  of  the  navigator  with  respect  to  the 
ellipsoid  and  the  values  from  equation  I.5»  we  obtain  an  expression 
for  the  radius  vector  at  any  time  as  shown  below. 
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t       v 


(lb    +  ho  )cos<po    cos\o    +  ?%HBt  + 


R  = 


^8    du  dv 


0        0 


t        v 
(No    +  ho  )coscpo    sia\o    +  to-  N„  t   +      /        /      j\9    du  dv 

0        0 


_t        v 


(No  (1-e3  )   +   ho)sin<po    +  i§Net  + 


Ynb  du  dv 


0   0 


t   v 
+  (MB  )t  +  /    /   Qx^xNb)  du  dv 
0   0 


?   1.6 


The  last  two  terms  on  the  right-hand  side  of  the  equation 
carry  point  B  in  a  circle  about  the  earth's  spin  axis,  and  can  b< 
accounted  for  by  a  simple  rotation  about  the  Y8  axis*  For  an 
incremental  time  step  of  At  we  may  write 


R8(-nAt)  = 


cos(nAt)  -sin(QAt)  0 
sin(flAt)  cosfoAt)  0 
0  0  1 


Dropping  the  subscript  Ne  as  being  superfluous  once  the 

motion  of  point  B  with  respect  to  inertial  space  is  accounted  for, 

— ♦   — •    — ♦ 
and  introducing  the  symbol  AR  =  RN  -  Rg  ,  we  can  rewrite  equation 

1.6  in  terms  of  incremental  time  steps,  relating  the  position  at 

time  t  +  At  to  the  position  at  time  t  as  the  equations  which 

follow  indicate. 


:(j+1)  =  I.C-QAt).4(j)  ♦  AR(j)  +  R(.}At  ♦ 
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or 


Y° 

< 

AY0 

Y°~ 

Y* 

=  Ra(-OAt) 

H 

+ 

AY1 

+ 

I* 

Y8 

(.i+1) 

.4 

m 

AY8 

m 

.*■- 

At  + 


LiJ  (j) 


+  £(At)8 


Y°,  .  ,v  +  Y0,., 
(j+D      (3) 


? 


(j+1) 


+  Y1 


(j) 


Y8  ,  .  .  .  +  Y8  ,  . , 

(0+1)    (a) 


>  1.7 


Since  we  cannot,  in  practice,  measure  A  (Y^  ,  Y  ,  Y8  ) ,  we  must 

substitue  components  of  R  which  we  derive  from  A8  and  which  con- 
tain the  effects  of  non-normal  gravitation.   Thus  we  use  quanti- 
ties obtained  from  the  values  sensed  by  the  accelerometers.   In 
this  study  we  shall  construct  the  sensed  values  by  taking  true 
acceleration,  adding  the  true  (hypothetical)  gravitation,  and  then 
subtracting  the  normal  gravitation.   Substituting  the  quantities 
so  derived  into  equations  1.7  we  obtain  the  set  of  equations  shown 
below. 
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Y° 
I1 
Ya 


=  Rg(-flAt) 


(j+1) 


"<" 

~ay°~ 

Y°~ 

* 

+ 

AY1 

+ 

i1 

.4- 

(j) 

AX* 

(d) 

Y3 

- 

At   + 


(j) 


£(At)s 


3r    +    (\   +  Ag)cos(cp  +   4.)  cos(/V  +   n/coscp   +  Qt)    + 

I     +    (y   +  Ag)cos(ip  +   ^)  sinCV  +   ti/coscp  +  nt)    + 
Y8    +    (y  +  Ag)sin(cp  +   O 

+  Qa  (N  +  h)cos(p  cos(/V  +  nt) 


+  08  (N  +  h)coa<p  a±n(/V  +  nt) 


•       •       e 


(J) 


£(At)! 


i°    +    (y  +  Ag)cos(q>  +   O    cos(\.  +   n/costp  +  fit)    + 
Y1    +    (y   +  Ag)cos(<p  +   O    sin(,V  +   n/coscp  ♦  At)   + 
_Y8    +    (y  +  Ag)sin(<p  +  5) 

+  Qa  (N  +  h)cas<p  cosCV  +  fit) 


+  Q8  (N  +  h)cos<p  sin(V  +  Qt) 


•        •        o 


J(j+1) 


>  1.8 


-  £(At)s 


-  £(At)s 


(VN    +Oa(NN  +  hN))coscpN  cos(,\*  +  Qt) 

Wm    +  n8(NN  +  hw))cos(pN  sin(,V*  +  CH) 
Yn  sincpn 

(Yn    +  n8  (Nn  +  hN))coa<pN  cos(\n  +  CH) 

(Yn    +  ^(Nn  +  hH))cos*pN  sin(\*  +  Qt) 
Yn  sin<pN 


(j) 


(j+1) 
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where  At  is  sufficiently  small  that  linear  interpolation  between 
the  components  of  gravitation  at  t/.v  and  t,  .   .  v  is  permissible. 

Xn  actual  navigation,  where  a  computer  is  allocated  to  the  navi- 
gation task  and  position  is  updated  several  times  a  second,  an 
approach  similar  to  the  above  is  acceptable;  however,  rather  than 
use  a  small  time  increment  and  compute  many  points  per  second 
which  are  not  needed  in  this  study,  the  effect  of  curvature  in 
the  quantities  to  be  integrated  will  be  considered <> 

For  the  purposes  of  this  study,  the  true  initial  velocity 
and  the  true  acceleration  of  the  navigator  with  respect  to  the 
earth  will  be  zero.   The  latitude  and  geographic  longitude  rates 
of  the  navigator's  position  will  be  the  result  of  variations  in  the 
gravitation  vector,  and  should  be  small-   Equation  1 08  can  then 
be  written,  for  use  with  a  larger  incremental  time  step  in  the 
form  shown  below « 


7k 


Y° 

% 

AY° 

I1 

=  Ra(-fiAt) 

n 

+ 

AY1 

Y8. 

(  "Ul   } 

A 

(  A  \ 

AY3 

(j) 


Y° 


At   + 


(j) 


+  ^(y  +  Ag) 


cos(<p   +   4)   Av.Val. 


cos(/V  +    ii/coacp   +  fit) 


cos(<p   +    4)   Av.Val.    sin (\.  +    ri/cosip   +  fit) 
sin(<p   +   £) 


(At)s 


+  rtl8  (N  +   h)cos<p 


Av.Val.  |cos(\.  +  fit) 
.  |sin(\   +  fit) 


Av.Val 
0 


(At)3    - 


Si. 9 


Av.VaX.  I^Cym  +  C5a  (Nn  +  hN))costpN  cosC\N  +  Ot) 
Av.Val.  |(Yn  +  0a(NN  +  hN))costpN  sinCA^  +  fit) 
Av.Val.      YNsinipN| 


(At)5 


The  most  complex  of  the  average  value  expressions,  (yn  + 


f}3  (NN    +  hN  )  )cosxpN 


C  OS  "1 

(Xn  +  fit)  can  be  written  as  (fj.  (<p,h)  + 

SID.  J 


fa  (<p,h))f3  (<p)f4  (\.  +  fit)  =  Fi  ((p,h)f4  (\.  +  fit).   Since  <p  and  h  are 
orthogonal  to  /\.  +  fit,  we  know  that  the  average  value  of  the 
product  Fx  (tp,h)fi  (\,  +  fit)  is  the  product  of  the  average  values. 
The  quantities  <pN  and  hN  will  change,  at  least  initially,  very 
slowly  because  of  the  initial  conditions  imposed.   We  may  take 


75 


advantage  of  this  fact  and  write  the  average  value  of  Fi  (cp,h)  as 


Av.Val.    (Yn    +  fJ3  (Nn    +   hN))cos<pN 


=  2^Vn    +  n*  ^Nn    +  h*  ^co 


scpN 


^[(Yn    +  ^(Nm    +   hN))cos<pN 


(o) 


(j+1) 


,    and 


Av.Val.    |^YNsin(pNj  =  ^(yn  sin<pN  )  ^    ^    +  ^(yn  sin<pM  )  (*+1) 

For  the  average  value  of  the  longitude  term  we  use  the  fact 
that   the  average  value  of  some   function,    f(x),    over  an  interval 
from  a  to   b  equals 


b-a 


f(x)   dx. 


A\.  +  C5  At 


\.  +  Qt  +  A\  +n  At 


\.  +  nt 


sinx  dx     = 


A/V  +f>  At 


-   cos(/V  +  fit  +  A\   +OAt)   +  cos(\.  +  fit) 


A\,  +0  At   L 


cos(X.  +  fit)    -   cos(/V  +  nt)cos(A\   +  QAt)    + 


+  sin(?u  +  nt)sin(AV  +  C3At) 
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l\.  +n  At 


cos(,V  +  CH)(1  -   cos(A\   +0  At) )    + 


+  sin(X.  +  f3t)sin(A\.  +n  At )      s 


=  cos(/V  +  nt) 


a\.  +  n  At     (a\.  +  n  At)3 

21  4l  + 


+   sin(\   +  CJt) 


1  - 


(AX.  +QAt)3        (A,V  -t-fl  At)* 


3i 


51 


-•••] 


Because   of   the   initial   conditions  imposed,   A\~0,    so 
A\.  +CJAt^rnAt=  0.005  for  At   =  60  seconds. 


(0At)3«r  10"7 


(ft  At)4-  10"9 


(QAt)3_    ..-8 

-TT--10 


(ftAt)< 
51 


10 


-11 


Since  the   error   of  truncation  of  a  series  of  this  type  is 
less  than  the   first  neglected  term,   we  are   safe   in  going  only  to 
the   cube  of  A\.  +0  At. 

Av.Val.feinU  +  Ot)]    =  sintV  +  Clt)(;j)[l  -   (A/V  ^"At)8  ]  + 


+  cos(\.  +  nt) 


(o) 


A\  +  QAt 


[*- 


(A\.  +  QAt)a 
12 


Similarly, 


Av.Val.TcosCX.  +  nt)J    =  cosCV  +  nt)(j)[1  "  — ^~Z 


(A\.  +  OAt)a 


•    f*        rt*\          A\   +C1  At 
-   8in(\.  +  fit)/  .v   2 


1  - 


(A,V  +  fiAt)a 


12 


] 
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Substituting  the  average  values  of  the  longitude  and  latitude 
functions  into  equation  1.9  we  obtain  the  final  equations  for 
rectangular  coordinates  of  position  in  an  inertial  reference  frame 
as  shown  below* 


(j+1) 


cosQftt      -einQ  At      0 


sin f}  At        cosfiAt      0 


0 

'4 

~AY°~ 

0 

4 

+ 

AY1 

1 

J*. 

(j) 

_AYa 

Y1 


At  + 


J  (d) 


(d) 


cos(*   ♦   O  (cosU   ♦  ^U  +  Qt)  [l  -  -&^-3 


sin(<p  +   £) 


-   sin 


a  +  -JL.  +  Qt)  &|±  [l  - 
cosu>  2      L 


(flAt)* 


coscp 


+  cos(\.  +    J*  '  +  nt 

COSU) 


)Qii[i- 


12 

(n  At)a 


12 


|(V  +  Ag)(At)a    ♦ 


+  if^CN  +   h)cos<p   (At)3 


cos 


a»flt)[i.lf]. 

sin(V  ♦.  0»>  [l  -  *£*■]    ♦ 


IolO 
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-   sinU  +  Ot)  '■=-—  II  - 


2 


♦   coa(V  ♦  Qt)  ^ 


1  - 


(nat)8 

12 
(ft  At)8 


12 


l-.Vi 


-  t(At)1 


([(Yh    +na(NN    +   hN))coscpNl  (    j    + 
([(Ym    +«a(NN    +   hN))cos<pMl  /  .j    + 


+    [(Yn    +08(Nm   +  *"^WB»ii]-(j*i)]  fC08^M   +  Qt)[l 
+    [(Yn    +Oa(NM    +  h«i))C08<Pii](j+1))  (sin(^    +Ot)[l 


•    •    • 


-   Uh    tatit)'}-  .inU,    *  Qt)  &i_|flAi  [i  . 

.  (ft,  *flAt>*j,  e„a,  +  nt)  ^_±n&t  ^  . 


•      •      • 


(a\m  +  nAt)8 

12 


(Aj^    +  QAt)8 
12 


APPENDIX  II 

DEVELOPMENT  OF  FORMULAS  FOR  SENSORS 

ROTATING  WITH  THE  EARTH 

In  the  case  where  the  accelerometers  maintain  a  fixed  orien- 
tation with  respect  to  the  earth,  the  X  reference  system  is  used, 
where  X  is  identical  to  Y  at  time  t  =  0.   Since  the  reference 
frame  is  rotating  we  must  consider  the  effects  of  this  rotation 
on  the  measured  accelerations.  The  accelerometers  continue  to 
sense  accelerations  in  the  inertial  frame,  resolved  into  compo- 
nents along  the  X  axes,  rather  than  measuring  accelerations 
relative  to  the  X  reference  frame. 


r  =  R  +  p 
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•  •  •» 


pT    =  nxp    +   pi 


rY    =  Rf    +  nxp        px 


=  *»    +   dt (Qxp  )y    +  "XP*    +  P« 

^(nxp)Y  =  nxp  +  nxpY 


=  nxp  +  nx(nxp)  +  nxp, 


••  •• 


rY  =  Ry  +  nxp  +  nx(nxp)  +2nxpx  +  px  =  a 

In  our  case,  where  the  X  and  T  systems  have  a  common  origin, 

and  where  the  rotation  rate  of  the  X  system  is  constant  at  n,  we 

may  write 

•  »• 

a  =  n*x(n*xp)  +  2n*xpx  +  px  n.i 

As  -in  Appendix  I,  we  may  write 

!».-■-=  A  -  H  II. 2 

Although  A  is  the  acceleration  in  inertial  space,  it  need 

not  be  sensed  along  the  7° ,  Y1  and  Y8  axes.   In  the  present  case 

we  will  sense  along  the  X° ,  X1  and  Xa  axes.   The  gravitation 

vector  can  be  expressed  in  this  system  as: 

GM  =  -  (ym  +  na  (Nm  +  ha  ) )  cos<pM  cos\*  i 

-  (Yn  +  Q3  (Nn  +  hN))  cos(pM  sii&N  j  II. 3 

-  Yn  sincpN  k 
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The  expression  for  actual  gravitation  will  be  written  in  the 
same  manner  as  was  done  in  Appendix  I,  taking  advantage  of  the 
fact  that,  for  this  study,  the  navigational  experiment  is  being 
performed  at  or  near  an  equipotential  surface;  therefore,  gravity 
anomalies  will  be  used  rather  than  gravity  disturbances. 

Q  =  -  Uy   +  Ag)  cos(tp  +  £)  cos(\  +  n/cos(p)  + 
+  fiaNcosq>  cos\.   i 

-  \(y   +  Ag)  cos((p  +  £)  sin(\  +  n/coscp)  +  ZI»4 

+  Q  Ncosip  sxn/V   j 

-  I  (Y  +  Ag)  sin((p  +  O  k 

Putting  equations  II. 1  and  II. 3  into  II. 2,  we  write  the 
equation  that  is  used  to  solve  the  navigation  problem  from  the 
sensed  acceleration: 

a9  =     [-  n»xS   -20&  +  xS  * 

+   (yn    +  A3  (Nn    +  hH))coa<pH    cos^hJ    i 

+  [-  naxji  +  2nift  +  %  + 

L  II. 5 

+    (yn    +  ^a  (NN    +  hN))cos<pM    sin**!    "j 

♦  [  Xs,  +  Y»  sin(pM  J  k 
To  obtain  the  hypothetical  expression  for  sensed  accelera- 
tion, we  write  the  same  basic  equation,  but  for  the  theoretical 
position,  velocity  and  acceleration  of  the  navigated  instrument, 
and  under  the  influence  of  hypothetical  gravitation.   To  obtain 
this  result,  we  put  equations  II. 1  and  II. k   into  II. 2  and  obtain  the 
results  shown  below. 
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As    =      [-  Q3XP    -  2QX1    +  X?    +  OaNcos<p   cos,V  + 

+    (y   +  Ag)cos(<p  +   O    cos(\.  +   Tj/costp)J    i  + 

+  f-  QaX1    +  2  OX?    +  X1    +  OaNcoscp  sin/V  + 
L  II. 6 

+    (*Y   +  Ag)cos(cp  +   O    sin(\.  +   ri/coscp)     "j   + 
+  [    X?    +    (y   +  Ag)sin(<p   +   £,)J   k 

•  M 

-*  — * 

Setting  px    =  px    =  0,   and   equating  terms  between  II. 5  and 
II. 6,   we   obtain 

X$    =  -  na(X°    -  XS)    +2QX]i    +  C^Ncoscp  cos\.  + 

+    (y  +  Ag)cos(q>  +   5=)    cosCV  +   ri/cos(p)    -  II. 7a 

-  (yn    +  fla(NN    +   hN))cos(pM    cos\* 
X}j    =  -  C^CX1    -  Xi)   -20X5    +  naNcos<p  sinV  + 

+   (Y  +  Ag)cos(cp  +  £)   sin(V  +  n/cos(p)   -  II. 7b 

-  (yn    +  08  (NH    +  hn))coBtpH    sin\w 

X*    s  (y  +  Ag)sin(tp  +   O  -  YN8in<pN                                                    II. 7c 

In  view   of   the  fact   that 

X9    *  Ncos(p  cos/V  X$    =   (NN    +   hN)coa<pN    cos\M 

X1    =  Ncostp  sinV  XJi    ■   (NN    +   hrJcoecpn    sinA* 

equations  II. 7  reduce  to  the  following  expression. 
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ss 

ii 

a 

=   20 

-i°« 

ii 

0 

+  (V  +  Ag)cos(tp  +  O 


cos(V  +  Tj/cOS(p) 
sin(,V  +  Ti/cos<p) 


tan(<p  +  4) 


II. 7' 


-  YnC08<Pn 


coa\„" 


sinV* 


tarnpN 


i; 


Although  the  determination  of  the  XN  components  involves 
more  computations  than  did  the  determination  of  the  components  of 
XN ,  the  integration  which  follows  is  still  straightforward,  and, 
since  the  gravitation  vector  for  a  point  remains  fixed  with 
respect  to  the  rotating  coordinate  axes  (accelerometer  axes),  a 
simple  mean  may  be  used  to  obtain  the  average  acceleration  over  a 
time  increment  in  place  of  the  more  elaborate  methods  used  in 
Appendix  I. 


p*  =  p*o  +  po  t 


t   v 

+    /    /   p*  du  dv 
0   0 


This  equation  may  be  written,  incrementally,  in  a  manner 
similar  to  that  used  in  Appendix  I,  as 


Sk 


X°N 

A. 


(j+D 


"xg" 

r&" 

xi 

+ 

i& 

.* 

(j) 

.xS. 

At  + 


(j) 


+  |nUt)3 


o 


+  -r(At)a(v   +  Ag)cos(<p  +   £) 


cosCv  +   T^/coscp)" 
sin(\    +   q/coatp) 
tan(<p  +   4) 


>  IIo8 


-  ^(At)avN^  .jCOS(pM^^ 


cos/V>j 

sin\N 
tarupN 


(j) 


-    ^t)aVN(;j+l)C08VN(;j+l) 


cos\N' 

sin\M 

tancpN 


(j+1) 


APPENDIX  III 

DEVELOPMENT  OF  FORMULAS  FOR  SENSORS  ORIENTED 

TO  THE  ELLIPSOIDAL  NORMAL 

In  the  case  where  the  accelerometers  maintain  an  orientation 
with  one  axis  parallel  to  the  ellipsoidal  normal,  one  axis  in  the 
horizon  plane  in  the  meridian,  and  the  third  axis  in  the  horizon 
plane  in  the  prime  vertical,  the  Z  coordinate  system  is  used  pri- 
marily; however,  since  the  accelerations  are  sensed  in  inertial 
space,  it  is  necessary  to  transform  repeatedly  between  the  Y 
system  and  the  Z  system.   The  transformation  equations  used  are: 

%    =  Ri(-<p)  Ra(A)  RY  =  MTRy 

"+cos(p   cosA        +coscp  sinA        +sinq> 


MT 


-sinA 


+cosA 


-simp  cosA   -sin<p  sinA   +coscp_ 

As  before,  As  =  R  -  G.   Let 
R  =  B  +  V 
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IP    = 


cosA 
sinA 
0 


COS(p 

cosA 

CO  Sep 

sinA 

sincp 

B     =     Ne  cosip 


V°    = 


V     =     N(l-ea)    +  h 

R     =   B°B  +  V°V 

•  •  *         •  • 

"5    =  IPb  +  UPb  +  ?>v  -i-  ?>v 

=    (H  xlP  )B  +  IP  B   +    (<33w  X^  )V   +   V°  V 

•  •  •  •  •  o  •• 

R      =    (JUg  xlP  )B  +    («Jb  XEP  )B  +    (u>B  xlP  )B  +  IP  B  +  ]£  B  + 
•  •  •         •  >• 

+  ((uvx"^°  )v  +  (ujvx^  )v  +  (u^x?3  )v  +  rv  +  v°v 

«  • 

=   (aoe  XEP  )B  +    (ui)v  x"^0  )V   +  ou^  x C«%  xEP  )B  +  u)f  x(<uv  XV0  )V  + 
+  "3d  B   +  V°  V   +   2(tuB  xlP  )B  +   2(u)v  XV0  )V 


III.l 


III. 2 


MT1P    = 


MTV°    = 


coscp 

0 

-sin<p 

l" 

0 

0 
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B  = 


Ne4 

2W3 


co sip  sin2(p  -  Nea    sirup 


(p 


Nea(l-ea)   sintp      • 
= J3 <P 


B  = 


Ne' 


y  e3(l-ea)sirup  sin2<p  -  Ne  ^"e    }    coscp 


Ne4(l-ea) 
W4 


sirup   sin2cp 


•a        Ne3  (l-ea  ) 

w3 8     v  v 


3Ne4     /,      a  n    •  o  Nea     ,n      3  x 

Ca     d-e*  )sxnip  sm2<p  -  —5-   (1-e"*  )cos(p 


(p     - 


Ne3(l-e3  )sin(p     •• 

¥ * 


•          Ne3  (1-e3  )      ,    ,      •         • 
=     2W3 sin2<P  «P   +  h 


V  = 


Ne4(l-e3)         .   a0  Ne4(l-ea)      .   8o  Ne3(l-e3)    3 

jjg5 sin82cp   +       ■   2Vtf4 sin32cp  +  —5 2  cos2<p 


<p3  + 


Nea(l-ea)  -     "       " 

+  -5 sin2cp  <p  +  h 


3Ne4(l-ea)      .   8o          Ne3  (l-ea  )    „      - 
* jj; sm82<p  +  jj, cos2<p 


"a 

<P        + 


Nea(l-ea)      .    _     ••       v 
+  -5 sxn2<p  <p  +   h 


W«     = 


0 
0 

n  +  i. 
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44 


U)y         X 


<p 

sinA 

-(J) 

cosA 

n 

+  &  . 

(p  sinA   +  (p  cosA    (Q   +  £.) 
•<p  cosA   +  <p  sinA    (ft   +  /V) 


MT(1PB  +  V°V)    = 


v        Nea(l-ea  )      .   a      .a 
W5 sin'cp  <pa 


Nea  (l-ea  )    sin2(p  3ea       ,   a  •8 


W3 


>    III. 3 


Nea(l-ea)      ,   a      « 
♦  $? sin  <p  (p 


uSexlP    = 
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cosA        sinA 


n  +  ,Y 


-sinA    (0   ♦  /V) 
cosA    (G   +  ,V) 
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UOy 


X  V°    = 


tp  sinA 


-<p  cosA  0   +  \ 


cosip  cosA        costp  sinA        simp 


*     - 

-  sincp  cosA  <p  -   cosip  sinA    (0  +  \.) 

-  sirup  sinA  <p  +  coscp  cosA   (0  +  \J 
cosip  q> 


MT  (tt,  xlP  )    = 


o  +  \. 


MT  (ujv  x"tf°  )    = 


cos<p  (n  +  \.) 


«p 


2MT 


(3fcxlP)B   +    (^,XV°)V]    = 


_  2  Nea(l-e    )    sin3  f   '(Q  +  Jj    +   2  cos<p  h(Q  +  *J 


Nea  (1-e2  )    sintp  COS(p  <Ja    +   2  h<p 
W3 


>III.Jt 
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tUeXET    =     0  0  \ 


cosA        sinA        0 


-  sinA  \ 
cosA  \ 
0 


m 

MT  (ujg  xl0  )    = 


0 
\ 


~*       -to 
UU^XV0     = 


"TF  "^  ^^ 

«• 
..  •  .  • .  ••  •  • . 

sinA  cp  +  cosA  cp(Q   +  \.)        -cosA  cp  +   einA  cp(Q   +  \J        \. 


coscp  cosA 


coscp  sinA 


simp 


-  sincp   cosA  <p  +   sincp  sinA  cpui   +  \J    -   coscp  sinA  \ 

-  sincp  sinA  cp  -   sincp  cosA  cp(Q   +  \)    +  coscp  cosA  \. 
+   coscp  <p 


MT(w,X?>)    = 


-   sincp  cptft   +  \)   *   coscp  \, 


<P 


mt  [(3bx1P)b  +  («,x?>)v]  = 


0 


(N  +   h)    coscp  'i.  -    (N(l-e3)    +   h)    sirup  cp(Q   +  £j 
(N(l-ea)    +  h)   'cp 


u>b  X  (<i  X IP  )    = 


0 


k 


-   sinA    (0   +  \)        cosA    (Q   +  \.)        0 
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>  III. 5 


-  cosA   (0   +  \Ja 

-  sinA    (ft   +  \Ja 
0 


MtC^xCIubxIP  ))   = 


-   coscp    (Q  +  \.)a 
0 
sincp   (Q   +  \.)a 


iu,  X  (u^  XT?0  )    = 


(p  sinA 

-   sincp   cosA  cp  - 
-   coscp   sinA    (ft   +  V) 


-  cp  cosA 


-   sincp  sinA  <p  + 
+  coscp  cosA    (Q   +  /V) 


k 

Q  +  i 

coscp  cp 
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(JO,  X  (<3jy  X$°    )      = 


-  cos<p   cosA  <pa  +   sincp  sinA  cp(Q  +  \.) 

"a  •  • 

-  coscp   sinA  (p  -   sincp   cosA  cp(Q  +  \) 

-  sincp  (pa 


•      _ 

-   coscp  cosA   (Q   +  \) 


-   coscp   sinA    (fl   +  \) 


^a 


•    •    » 


MT  (^xduyXV0))   = 


-  cp8    -  cosPcp   (n  +  \.)a 

-  sirup  cp(C2  +  \) 

_  +   sincp   coscp    (CI  +  \.)a 


MT    (Ub  x  (3b  x IP  )  B  +  wv  x  (<*>»  X V°  )vl    = 


-  (N(l-ea)    +   h)   <pa    -    (N  +   h)cosacp   (Q   +  ,Y)a 

-  (N(l-ea)    +  h)    sincp  <p(fl   +  \.) 
+    (N  +  h)    sincp   coscp    (Q   +  \.)a 

r  — ♦ 

Equations  III. 3  through  III.o  provide  the  components  of  R 

•• 
that  go  into  equation  III. 2.   To  determine  R  as  a  function  of  the 

sensed  acceleration,  it  is  necessary  to  obtain  the  components  of 

the  gravitation  vector  also.   From  considerations  similar  to  those 

of  Appendix  I,  we  will  derive  the  gravitation  vector  from  the 
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gravity  vector,  since  we  are  considering  navigation  on  or  near  the 
reference  ellipsoid  in  the  ideal  case. 

-1  +  nx(Qx(N  +  h)l°)  =  -  v 


MT  (-G)    = 


Y   +    (N  +   h)cosa<p  Cls 


-    (N  +   h)    sin<p   coscp  C3a 


III. 7 


A,    = 


s& 


3 


e     rn(i-e° ) 


H 


w* 


+   h 


(p8    -    (N  +   h)cosa<p  \(20  +  \.)    +  y 


(N  +  h)costp  £  -   2[        ~6        +  h|    sin<p 


'] 


a  \    jb 


fN(l-ea)        .I  -        3N(l-ea  )    ea  „      -a  •"• 

[— ^5 +   hj   cp  +  ^-^g Jr  sm2<p  <pa    +   2hq> 


+  2  coscp  h(Q  +  \.) 


+  (N  +  h)  simp  coscp  \(2Q  +  ,\.) 


N(l-ea  ) 

In  view  of  the  fact  that  s is  the  radius  of  curvature 

VT 

of  the  ellipsoid  in  the  meridian,  M,  we  may  rewrite  the  equations 

for  the  sensed  accelerations  in  the  more  compact  form  shown  below, 


s* 


IS 


2ft 


"a 


#  •  g 

h  -    (M  +  h)cp3    -   (H  +  h)cos8<p  \.(20  +  \)   + 
(N  +   h)ct5S(p  \.  -   2(M  +   h)sin(p  <p(fl   +  £)    + 


/  \**  Me  •»  °  • 

CM  +   h)<p   +   3  ■—=-   sintp   coscp   <p      +    2hcp   + 

w 


+   Y 

+  2cos<J>  h(n   +  \.) 

+    (N  +   h)    sincp   coscp  \(2Ci  +  \.) 


III. 8 


The  above  equations,  when  evaluated  at  the  navigator's 
position,  permit  the  determination  of  h,  \.   and  <p  from  the 
sensed  accelerations,  under  the  assumption  of  normal  gravity. 
These  equations  are  equivalent  to  the  kinematic  equations  of 
motion  presented  in  standard  works  on  inertial  navigation,  but 
these  equations  are  in  terms  of  fundamental  parameters  used  by 
geodesists  (M,  N,  ea  ,  tp,  \.,  h)  rather  than  in  terms  which  have 
more  general  applicability,  and  which  are  used  to  fit  any  situ- 
ation (radius,  geocentric  latitude,  longitude),  but  which  must 
then  be  transformed  to  obtain  the  quantities  which  the  navigator 
desires. 

As  in  Appendixes  I  and  II,  a  problem  will  be  considered  in 
which  a  navigation  system  is  used  which  is  initially  in  perfect 
alignment  with  the  coordinate  system  at  the  geodetic  position  of 
the  observer,  but  in  which  the  gravity  vector  does  not  coincide 
with  normal  gravity.   The  centrifugal  effects  given  in  equation 
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III. 7  will  remain  unchanged,  since  they  are  a  function  of  position 
only,  but  the  net  gravity  effect  will  change  to 


(g)z    =  Ri  (-<p)Ra(A)R3(-  A  -  n/coscp)Ri  (<p  +  £) 


V  +  Ag 


Using  this  relationship,  and  tacitly  assuming  that  we  have  no 
geoidal  undulation  with  which  to  contend,  it  is  possible  to 
rewrite  equation  III .8  for  the  sensed  accelerations  in  the  hypo- 
thetically  true  case  as 


"*. 


••a 


h  -    (M  +  h)ia    -    (N  +  h)cos8(p  ,V(20  +  \) 

»•  •  • 

(N  +  h)coe<p  X.  -   2(M  +  h)sintp  <p(0   +  \)   + 


Me 


(M  +   h)(p  +  3  -38"  sin«p  cos<p  tpa    +  2hq> 


•    •    • 


+  2cos<p  h(CJ   +  \) 

+   (N  +  h)sin<p  coscp  \(2Q  +  \J 


^111.9 


+    (v   +  Ag) 


sintp  sin((p  +   4=)    + 


•    •    • 


coscp   smv.<p 


(«  +   £)- 
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+   coscp  cos(<p   +   5)    cosA   cos(A   +    ri/cosq?)    + 

-  cos(cp  +   ^)    sinA   cos(A   +   ti/coscp)   + 

-  sincp  cos(<p  +   O    cosA   cos(A   +   Tj/cosip)   - 


+   cos<p   cos((p  +   4.)    sinA   sin(A   +   Tj/cos<p) 

+   cos(<p   +   O    cosA   sin(A   +   r\/cos(p) 

-   sincp   cosCcp   +  O    sinA   sin(A   +   n/coscp) 


J 


In  the  cases  which  have  been  considered  in  the  first  two 
appendixes,  the  navigator  has  been  stationary  with  respect  to 
the  earth,  so  that  h  =  <p  =  \.   =  h  ■  q>  =  \.  s  0.   In  this  situation 
equation  III. 9  simplifies  to 


zl 


z? 


a     (Y    +    Ag)cOs((p    +    O    MT 


cos(A  +  n/cos<p) 
sin(A  +  q/coscp) 


tan(<p  +  O 


III. 10 


If  the  stable  platform  were  to  maintain  its  orientation  with 
respect  to  the  normal  to  the  ellipsoid,  and  to  rotate  such  that 
the  Z8  axis  remained  aligned  in  the  meridian,  the  above  quantity 
would  remain  constant.   Unfortunately,  the  rate  of  driving  the 
platform  about  the  various  axes  is  a  function  both  of  indicated 
velocity  and  indicated  latitude,  and  the  axes  themselves  are  con- 
stantly moving,  making  the  problem  difficult.   The  rotation  rate 
about  the  instantaneous  Zx  axis,  which,  for  the  navigator,  is 
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defined  as  east  regardless  of  its  actual  orientation,  is  cpN .   The 
rate  about  the  instantaneous  Z  axis  is  (ft  +  /\n)cos(pN,  and  that 
about  the  instantaneous  Zi  axis  is  (fi  +  /X^Jsin^^,.   Since  these 
rotations  are  not  vector  quantities,  the  order  in  which  each 
rotation  is  accomplished  is  important;  however,  these  rotations 
are  accomplished  simultaneously  in  practice  and  in  infinitesimal 
increments. 

Let  us  consider  that  during  one  time  increment,  dt,  the 
transformation  matrix,  MT  ,  is  modified  by  the  transformation 

My  ,         v  =  R0((n  +  \^)sin<pN  dt)Ri(-<pN  dt)Ra((Q  +  \H  )cos<pM  dt)MT^.^ 

This  expression  would  become  quite  unwieldy  when  expanded;  how- 
ever, if  we  limit  dt  such  that  we  may  consider  the  cosines  of  the 
resultant  angles  to  be  one,  the  sines  as  the  angles  in  radians, 
and  the  product  of  two  angles  to  be  negligible,  the  order  of 
multiplication  becomes  unimportant,  and  the  relationship  may  be 
approximated  by 

MT(j+1)    =  MeMT(j) 


where 


^    = 


+1 


+ANcos<pN    dt        (pN    dt 
-ANcoscpN    dt        +1 


« 

-cpN    dt 


AN  sincpN    dt 
-AH  sintpN    dt        1 


where  AN    =  Q   +  in »      At  any  time   t  we  may  write 
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ZS 


•El 

z1. 


i% 


t 

w 

v=0 


=    (V  +  Ag)    TT    M6(v)MT(0)cos((P  +   V 


cos(A  +   Tj/coscp) 
sin(A  +   ij/coscp) 


tan(cp  +  £.) 


III. 11 


The  evaluation  of  equation  III. 11  can  be  extremely  tedious 
because  of  the  small  time  step  required  in  the  evaluation  of  the 
matrix  Mg(v).  We  may  rewrite  that  equation  in  terms  of  individual 
rotation  matrices  in  the  vollowing  manner: 

%9  ■  ■  TT   Ri(-<PNidti)Ro((n  +  \Ni)sin(pNidti)Ra((n  + 

+  /Ui)cos(pM;Ldti)   Ri(-(p)Ra(/V)  x  III. 12 


ff  ("HtC-Qdt,) 


gx 


where  the  product  Ra(-Qdt.)  transforms  the  gravity  vector  from 
X  space  to  inertial  (Y)  space.   The  product  Ri(-<p)R»(\)  trans- 
forms it  to  the  Z  space  as  defined  at  time  t  =  0,  and  the  long 
product  which  makes  up  the  remainder  of  the  expression  transforms 
to  Z  space  as  defined  at  time  t,  which  corresponds  to  the  instan- 
taneous sensor  axes. 

The  two  rotation  matrices 


Ro(fo  +  £*)sin<pN  dt)R8((n  +  \*)cos(pM  dt) 


represent  nothing  but  the  component  rotations  of  (0  +  \#  )  about 
the  vertical  and  north  axes.   These  can  be  replaced  by  a  matrix 
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Rs((n  +  \*)dt)  placed  where  the  axis  of  rotation  is  either  Xa  or 
Y3. 


C  J5  (0  +  in  ) 


(0  +  /Vm  )coscpN 


(fi  +  \* )  sin(pN 


Thus,  equation  III. 12  may  be  written 


%*    =  TJ   ■»(^1it1)J  Ri(-<p)R3U)  x 


*  IT 


RaUfl  +  in  Jdt.)R8(-0dt.)l  g 

J    J        J  J   x 


=  TT  [H1(-JMidti)]  Ri(-(p)R8(\J  Tp  [HgCi^^dt^)]  £x     III. 13 


Since  the  product  of  rotation  matrices  about  the  same  axis  is 
merely  the  rotation  matrix  for  the  sum  of  the  individual  argu- 
ments, equation  III. 13  may  be  written  as 


Rzs  =  Ri 


-<pNdt 


L  0 


Ri(-<p)Ra(V.)Ri 


\*dt 


L  0 


St 


=  Ri(-<PN  )Ra(?c«  )(y  +  Ag)cos((p  +  O 


COs(\.  +  T]/cOStp) 

sin(\.  +  Tj/cos(p) 
_tan(<p  +  4.) 


III. Ik 
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Combining  equations  III. 8  and  111.14,  we  obtain  expressions 
for  the  observed  accelerations  of  the  navigator  in  height,  long- 
itude and  latitude  as  follows. 


hN 


(N  +  h)NcostpN    \N 
(M  +  h)NtpN 


+    (M  +   h)N    fa    + 

+    2(M  +   h)Nsin<pN    cpN  (Cl   +   \N  ) 

-   3MNeasin2cpN    fa/ 20%    - 


+    (N  +  h)NcosacpN    \*(2n  +  Xh  )   -  yn 

-  2cos<pN    hN  (C3   +  \u  ) 

-  2hN<pN    -    (N  +   h)NsincpN    cos<pN   \*  (20  +  \*  ) 


+  Ri(-<pN  )Ra(X«  )(y  +  Ag)cos(cp  +  £j 


cos(\.  +   Tj/costp) 
sin(/V  +  n/coscp) 


tan(cp   +  5J 


> III. 15 


**    av   =   2   ^  +  Ag)cos(<p   +   4)   x 


-   sin\N    cos(/V  +   Ti/cos<p)    +   cos\N    sin(\    +   ri/cosq?) 
(NN    +    hw  )coscpN 


(j) 


-   sin\fl    cos(\i,  +  n/cosy)    +   cos\w    sin(\   +   m/cosq?) 


+  ' 


(NN    +   hN  )cos<pN 


(j+1) 
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(MN    +   hN  )sin(pN    q?N(n  +  j*  )    -   costpNhN  (fl   +  \^  ) 
(NN    +   hN  )coscpN 


(j) 


KMfj    +   hn)sincpN    cpN(Q  +  \H  )    -   cos(pN    hN  (fl   +  \^  )  ' 
L  (NN    +   hN  )coscpN 


(j+1) 


Collecting  terms,  and  putting  hN  ■  h*,  =  0,  we  obtain 


In 


av 


(y   +  Ag)cos(cp  +   EJ 
2 


sinC/Vw    -  \   -   Ti/cos(p) 

NNCOS(pN 


(o) 


rsin(\w    -  \   -   n/cos<p)1  +  |  MM  sin<pN (pN  (0   ♦  A^i 

L  NNcoscpM  J/  .  ^t  \     '  L  NNcos<pN 


(j+l)J 


>  III. 16 


(j) 


MM  sin(pN    u)M  CO   +  \m  ) 

NN  COS(Pn 


(J+D 


<Pn 


v^g 


=   >W»     ■  7^n I"-   sin<pN    cos/\^    cos(<p  +   4)    cosCV  +   n/coscp)- 

2V.MN    +  fan  ;r  H^    L 

-  sintpM    sin\,j,    cos(q>  +   £)    sin(\.  +  Ti/cosq>)  + 


av  '  ( -i ) 


+   cos(pN    sin(<p  +   4) 


(d)    + 


2(m,  ;  Q^j 


-  sin(pN    cos/Vm    cosCcp   +   £)    cosCA,.  +  Tj/coscp)- 

-  sincpN    sin\w    cos((p  +   £)    sinCV  +   i]/cosq>)  + 
+  coetpM    sin(<p   +  ^^1(^+1)   " 


-  4(mh  !e3h.)(.)  [M»8in2*»  **•](»* 


hn(PM 


LMN    +   h„ 


(j) 
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[1\VL  +\)  sin2<PN  ^2n  +  £■)! 

(j) 


3ea 

MM*      ♦     hN) 


(j+tf 


MNSin2<pN    <p?/WS 


(j+D 


hn  <Pn     1 
IM„    +   hM  J 


(3+3) 


-U(M//hN)  sin2«*  ir<«0+  «•>] 


<Pn 


av 


Collecting  terms,    and  putting  hN    =  hN    =  0,   we   have 


=  2M ~   sin^N    cos(<p   +   4)    cos(\^    -  ,V  -   rj/costp)    + 


+   cos(f>M    sin(<p   +   4) 


(3)   + 


-   sinq>N    cos((p  +   O    cosCiif    -  /V  -   T\/costp)    + 


+   coscpM    sin(<p  +   ^b) 


J (5+1) 


>III.17 


-  Tp  [sin2(pM    <p*  /»n    (  j  )   -  [73^  sin2<pM    \*  (  2  fl  +  **  ^ .  ( j  )    " 

3e8   f  "a      a  1  T  N  *  • 


Once  the  average  accelerations  in  latitude  and  longitude  are 
computed,  the  determination  of  position  follows  as 


*(d+l)  =  9(j)  +  ^j)At  +  iUt)°  *av 


The  development  in  the  appendix  through  equation  III. 15 
would  indicate  that  computation  for  height  would  be  possible  in 
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addition  to  the  computation  for  latitude  and  longitude.   The 
development  was  not  completed  since,  as  the  results  of  the  two 
previous  sensor  configurations  show,  a  solution  for  height  in  an 
imperfectly  known  gravitational  field  leads  to  ridiculous  results. 
One  of  the  primary  objects  of  the  present  orientation  of  the  sen- 
sors is  to  isolate  the  latitude  and  longitude  sensors  from  the 
height  sensor,  recognizing  that  it  is  necessary  to  supplement  the 
inertial  data  with  information  from  another  source  for  height 
information. 


APPENDIX  IV 

DEVELOPMENT  OF  FORMULAS  FOR  NAVIGATION  ON  A 

HYPOTHETICALLY  UNDULATING  SURFACE 

In  the  first  three  appendixes  we  have  considered  the  effects 
of  anomalous  gravity  on  a  navigator  which  was  stationary  with 
respect  to  the  surface  of  the  eartho   Referring  to  equation  111*9 
we  see  that  by  having  the  navigator  stationary  we  have  eliminated 
a  large  part  of  the  equation,  which  could,  in  fact,  provide  a 
significant  contribution  to  the  sensed  value  of  acceleration*   In 
this  appendix  we  shall  consider  the  more  realistic  case  of  a 
vessel  moving  over  the  surface  of  the  earth,  wherein  it  is  sub- 
ject to  a  varying  effect  of  gravity,  and  where  all  the  factors 
which  contribute  to  the  sensed  value  of  acceleration  must  be  con- 
sidered* 

Let  us  consider  the  case  where  the  geoidal  undulations  may 
be  expressed  by 

<  ■  <°  [•"  %hr&  -  i]  [oos  ^f^  -  i] 

where  \o  »  \i  »  <(©  »  and  <pi  are  the  boundaries  of  the  tessera  in 
which  undulations  are  to  be  oomputed  and  ^Co  is  the  value  of  the 
undulation  at  the  center  of  the  area*   Using  the  relationships 

1  dC  1   *£ 

^="Rol?      and      n="  RcTsl?  dfc 
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we  may  write  analytical  expressions  for  the  components  of  the 
deflection  of  the  vertical  in  the  following  form. 

r        £a.  l\.„„  2n(\.  -  \o  )    ,1    2n        2n((p  -  <&  )  Ttr  w 

Rcostp    L  vtpi    -  cpo  ;  J\i    -  \o  C\i    -  \o  ) 

Although  4  and  tj  are  directly  dependent  on  the  variations  of 
the  undulation,  £,  there  is  no  such  rigid  dependency  between  £ 
and  the  gravity  anomalies o  We  may  start  from  the  principal  equa*- 
tion  of  gravimetric  geodesy 

and  see  that  since  T,  the  disturbing  potential?  has  not  been 
specified  in  a  form  with  a  unique  derivative,  we  are  relatively 
free  in  our  choice  of  the  hypothetical  distribution  of  anomalies. 
The  tests  described  in  Chapter  II  seem  to  indicate  that  the  mag- 
nitude of  gravity  is  relatively  unimportant,  so  long  as  we  have 
an  independent  means  of  measuring  height o   As  a  check  on  the 
importance  of  Ag,  we  shall  let  it  vary  in  two  different  ways, 
then  compare  the  resultso   In  one  case  we  shall  assume  that  Ag 
in  milligals  equals  Q   in  decimeters,  while  in  the  other  case  we 
shall  assume  the  negative  of  this  relationship  to  hold* 

Let  us  consider  the  instant  of  time  at  which  we  are  steaming 
at  a  speed  V,  and  at  an  azimuth  <*»   The  relationships  which 
relate  the  rates  of  change  of  latitude,  longitude  and  elevation 
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to  velocity  and  azimuth  on  the  surface  are 


a         _a  s% 


V  cos(£   +  Tl  )   COSOf 
♦  =     M  +  H  +  C " 


v   V  cos(4a  +  Ji-  sint» 
(N  +  I  +  C)cos(p 

h  =  V  sinUa  +  r?)% 


Since,  for  this  problem  in  navigation,  M,  N,  4,,  t|  ,  C,  tp,  V, 
H  and  Of  are  functions  of  time,  the  analytic  expressions  for  (p,  \. 
and  h,  which  are  needed,  in  equation  IIIo9,  get  quite  lengthy..   In 
order  to  avoid  the  use  of  a  very  small  time  increment,  we  will 
establish  the  simplifying  condition  that  V,  H  and  a   are  constant; 

o 

however,  the  fact  that  H  is  constant  does  not  imply  that  h  is 

o 

zero,  since  h  must  take  care  of  the  variations  in  ^°   Since  the 
establishment  of  analytic  functions  for  <p,  \   and  h  is  a  formidable 
task,  even  with  the  simplifying  conditions,  we  will  use  numerical 
methods  to  arrive  at  these  quantities <> 

o  o  o 

We  know  that  <p,  \.   and  h  are  well-behaved,  slowly  changing 
functions..   Using  standard  prediction  techniques,  we  will  fit  a 
curve  to  the  functions  at  the  four  most  recent  time  points,  then 
use  this  curve  to  predict  the  value  of  the  function  at  the  time 
point  next  to  be  considered,  and  the  point  immediately  thereafter 
The  difference  in  the  function  at  t,.  -v  and  *(j.i)i  divided  by 

the  time  increment  between  the  two  points  will  give  the  value, 
to  an  acceptable  order  of  accuracy,  of  the  time  derivative  of  the 
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function  at  t/.x.   The  imposition  of  an  iterative  cycle  if  rigid 
limits  on  the  difference  between  predicted  and  computed  positions 
and  functions  are  not  met  insures  that  we  will  not  use  values 
which  are  significantly  in  error  in  any  of  our  equations » 

As  a  set  of  equations  equivalent  to  equations  IIIolO  and 
III. 8,  we  may  write 


« 


i\ 


=   (v   +  Ag)cos(q>  +   2p   MT 


cos(Au.  +    Tj/coscp   +  Qt) 
sin(\>.  +   rj/coscp   +  fit) 


tan(<p   +   O 


h  -    (M  +   H  +  C)<J>3    -    (N  +  H  +  CJcos8^  i.(2Ci  +  \) 
(N   +   H  +  Ocosip  \.  -    2(M  +   H  +   C)sincp  <p(Q   +   \d    + 
(M  +   H  +   £)ip   +   3Me   sirup   costp  cp  /Vr    +   2hcp   + 


V        ■        © 


+  2cosip   h(0   +  \.) 

+   (N  +   H  +  C)sinip   C08<p  /V(20  +  \J 

hN  -  (Mm  +  hw  )(p8  -  (Nn  +  ht»)cos  ipN  k*(20  +  tai )  + 
(NN  +  h^cosipN  X*  -  2(MN  +  hN)sin<pN  <pn  (Q  +  kn )  + 
(MN    +  hN  )tpN    +  3MN  ea  sintpN    coscpN    ipn/Wn    +  2hmpn    + 


>     IV.  k 
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+  Yn 

o  o 

+   2cos(pN    hN  (Q   +  \tt  ) 

+    (NN    +   hM)sin<pN    cos(pN    \n  (20  +  \n  ) 

In  equation  IV. k%    hN  is  measured  from  the  geoid,  but  the 
navigator,  in  applying  it,  assumes  that  it  is  measured  from  the 
reference  ellipsoid  and  assumes  that  the  geoid  coincides  with  the 
ellipsoid.   After  going  through  a  derivation  similar  to  that  pre- 
sented in  Appendix  III  we  arrive  at  an  equation  which  is  the 
equivalent  of  equation  III.15<> 


*V 


(NN    +   hN  )cos(pN    \n 


.(MN    +   hN  )<pN 


h   -    (M  +   H   +  CV    ~ 

M 

(N  +  H  +  C  )cos<p  /V  - 

(M   +   H   +  C)°<P   +   3Measin2cp  <J>8  /2W3    + 


-  (N  +   H   +  C)cos3(p  \(20+  id 

-  2(M   +   H  +  C)sin(p  (p(Q    +  X.)    +   2cosq>   h(fi   +  ,\.) 
+   2h<p  +    (N  +   H  +  £)sin<p   coscp  \.(2Q  +  \J 


>      IV. 5 


+  Ri(-q>N  )RsUn  )(*Y   +  Ag)cos(cp  +   O 


cos(\.  +   T\/coscp) 


sin(\>.  +   Ti/coscp) 


_tan(q>   +   O 


(Ma  +  hN  )cpJ5  +  (NN  +  hN)cosa<pN  ^(20  +  /\u»  )  -  yH 
2(Mn  +  hN)sin(pN  <J»N  (n  +  \m)  -  2costpM  hM  (fi  +  \*  ) 
-3M*  e8  sin2<pN    <J>n/2Wn    -   2hN<pN    - 


6       e       a 


•       •       © 


©  o 

-    (NN    +  hN)sin<pN    cos(pN   \n(2G  +  ta  ) 
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Expanding  this  equation  we  get 


x  (v  +  Ag)cos((p  +  EJ   .  /,     /       ,x 

(N  +  H  +  C)cos(p  \.  -  2(M  +  H  +  C)sincp  <p(fl   +  j.) 
(NN    +  hN  )cos(pN 


•  •  •  o 

2cos(p  h(0   +  \.)    +   2(Mm    -t-  hN)sin^N    q?N  (Q   +  ;u  ) 
(NN    +   nn  )cos(pM 


2cos<om    hN  (fl   +  ^g  ) 
(NN    +   hN)cos(pN 


Collecting  terms,    averaging  over  a  time  increment,   and 
setting  H  =   hN    =  hN    =0,   we  obtain  the   final   expression  for   the 
average  longitudinal  acceleration. 
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'Vm 


av 


NN  coscpN 


(j) 


[g^V   +  Ag)cos(<p  +   £)    sinCY 


1  •• 

+    tl/cos<p   -  \*  )    +  -(N  +  C)cos<p  \    _ 

-    (M  +  C)sin<p  tp(n   +  \J    +   costp  h(Q   +  \)    + 
+  M,,  sin<pN    <pN  (ft   +  \M  )  J  ,    .      + 


~(Y   +  Ag)cos(<p  +   £)    sin(\, 
(jh-1)L 


NN  coscpN 


1  ** 

+   Ti/cos(p  -  \N  )    +  -(N  +  C)cos<p  ,V  - 

-    (M  +  C)sincp  <p(fl   +  \)   ♦   cos(p  h(fl   +  i.)    + 

♦  MNsin<pN    (J^CCJ   +  ^^(j+i) 


>     IV. 6 


•  •  r 

(M*    +   hM  )(pM    x   (y  +  Ag)cos(<p   +   5)   costpN    tan(cp  +   O   - 
-   tan<pN    cos(\   +   Tj/cos(p  -  \N  ) 
+  Isinacpbe'M^/W3    +   (N  +   H  +  C)i(2Q+  {)] 
+    (M  +   H  +  C)<P   ♦   2h<p  -   2hN<pN    - 


^sin2(pN 


3eaMN<JS/WJ!  +  (NM  +  hN)\„(2n+  \*  )] 


Collecting  terms,  averaging  over  a  time  increment,  and 
setting  H  =  hN  =  h*  =  0,  we  obtain  the  final  expression  for  the 
average  latitudinal  acceleration., 


Ill 


<Pn 


av 


2X» 


J(j) 


(V   +  Ag)cos(<p   +   4)    coscpN  |tan(cp   +  4)    - 


-   tantpN    cos(<Y  +   Tj/coscp  -  \*  )      +  2h<p  + 
+  Sif^MeS3/*8    ♦    (N  ♦  C)   V-(2n+  £)] 


+  NN?u(^n+  \N) 


(j) 


}   IVo7 


2MK 


(j+1) 


(Y   +  Ag)cos(<p  +  4.)   cos(pN  J  tan(tp   +  4.)   - 


~   tan<pn    cos(\.  +   q/coscp  -  \H  ) 


+   2h<p  + 


+  cm  +  o;  -  *^-[3%^a/*  ♦ 


nn  \*(2n  +  i«)] 


(d+i) 


Once  the  average  accelerations  in  latitude  and  longitude  are 
computed  by  means  of  the  above  equations,  the  determination  of 
position  at  any  time  is  the  same  as  that  developed  in  Appendix 
III. 
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